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FUSION RULES FOR THE LOGARITHMIC N — 1 SUPERCONFORMAL 
MINIMAL MODELS I: THE NEVEU-SCHWARZ SECTOR 

MICHAEL CANAGASABEY, J0RGEN RASMUSSEN, AND DAVID RIDOUT 


Abstract. It is now well known that non-local observables in critical statistical lattice models, polymers and 
percolation for example, may be modelled in the continuum scaling limit by logarithmic conformal field theories. 
Fusion rules for such theories, sometimes referred to as logarithmic minimal models, have been intensively 
studied over the last ten years in order to explore the representation-theoretic structures relevant to non¬ 
local observables. Motivated by recent lattice conjectures, this work studies the fusion rules of the N = 1 
supersymmetric analogues of these logarithmic minimal models in the Neveu-Schwarz sector. Fusion rules 
involving Ramond representations will be addressed in a sequel. 


1. Introduction 

Superconformal algebras have a long history in mathematical physics, being intertwined with the development 
of string theory through their role as infinitesimal symmetries of superstrings (see [1], for example), as well as 
appearing as extended symmetries of the scaling limits of certain lattice models [2]. These infinite-dimensional 
Lie superalgebras each contain a Virasoro subalgebra that quantises the infinitesimal conformal symmetries of 
the plane. Field-theoretically, the superconformal algebras are customarily parametrised by the number N of 
fermionic partners of the energy-momentum tensor. The simplest examples, after the Virasoro algebra itself 
(corresponding to N = 0) are the N = 1 superconformal algebras: the Neveu-Schwarz algebra [3] and the 
Ramond algebra [4]. 

After this debut in superstring theory and statistical mechanics, the structure theory for Verma modules 
over the Neveu-Schwarz algebra was quickly settled. A determinant formula for their invariant bilinear forms, 
originally conjectured in [5], was proven in [6] and the possible submodule structures were elucidated in [7] 
with the result essentially repeating that of the Virasoro algebra. In particular, all non-trivial submodules are 
generated by singular vectors, the dimension of the space of singular vectors in any weight space is at most one, 
and non-trivial homomorphisms between Verma modules are necessarily injective. 

By way of contrast, the structures of the Ramond Verma modules can be more intricate [8]: submodules not 
generated by singular vectors can exist, there can be up to two linearly independent singular vectors of each 
parity in a given weight space, and there can exist non-trivial non-injective homomorphisms. These features of 
the Ramond case are now relatively well understood, but serve as a useful toy model for the more pronounced 
difficulties that one encounters when investigating the Verma modules of the N > 1 superconformal algebras. 
Despite the fact that the latter algebras have fundamental applications to mirror symmetry [9], the AdS/CFT 
correspondence [10] and Mathieu moonshine [11], it is fair to say that their representation theories remain 
poorly understood at best. 

Our purpose with this article is to explore some of the representation-theoretic aspects of the Neveu-Schwarz 
algebra that pertain to logarithmic conformal field theory. The more challenging exploration of the represen¬ 
tations of the Ramond algebra will be addressed in a sequel [12]. Here, the qualifier “logarithmic” means that 
the underlying vertex operator (super)algebra admits modules upon which the Virasoro zero mode Lq acts non- 
semisimply. A collection of reviews on this topic may be found in [13]. The results presented here are motivated 
by a recent lattice-theoretic study, reported in [14], although logarithmic conformal field theories with super- 
symmetry have been discussed in the past, see [15-19] for example. This study proposes a conjecture, based on 
numerical evidence, that the scaling limits of a certain collection of integrable lattice models are described by 
logarithmic analogues of the N = 1 superconformal minimal models. Here, these logarithmic superconformal 
models are explored directly in the continuum, as conformal field theories. 
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In particular, we study the fusion rules of a certain collection of Neveu-Schwarz modules that we will refer to 
as (Neveu-Schwarz) Kac modules. This usage follows the nomenclature established in [20-22] for Kac modules 
over the Virasoro algebra. The physical relevance of these modules is that they are believed [14,23] to identify 
the scaling limits corresponding to an accessible class of integrable boundary conditions for the underlying 
lattice models. Because of this, the continuum fusion rules for Virasoro Kac modules have been well studied 
both via lattice approximations [24-27] and direct calculations [28-33]. In contrast, the fusion rules of the N = 1 
Kac modules only seem to have been studied from the lattice point of view [14], although those of the simple 
(irreducible) modules appearing in the N = 1 minimal models have received continuum treatments, see [34-37] 
for instance. 

We begin, in Section 2, by recalling the N = 1 Neveu-Schwarz algebra and reviewing the necessary structure 
theory of its liighest-weight modules and Fock spaces. The latter are the superconformal analogues of the 
well known Feigin-Fuchs modules that appear in the Coulomb gas free field realisation of the Virasoro minimal 
models. Their structure theory is detailed in [38]. This theory is important for the present investigation because 
we define the N = 1 Neveu-Schwarz Kac modules as certain submodules of Fock spaces, following [20-22] in 
the Virasoro case. 

When a conformal field theory is rational, meaning that the relevant modules of the underlying vertex 
operator algebra are semisimple and that only finitely many simple modules appear, the fusion rules may 
be efficiently computed using the Verlinde formula [39]. While rigorous proofs seem to require rationality, 
see [40], a lesson learned from more physical proofs [41] is that this formula follows from the deeper consistency 
requirements of conformal field theory and so should be valid, in some form, more generally. In Section 3, we 
utilise a continuous version of the Verlinde formula to compute the characters of the fusion products of the Kac 
modules. This approach was pioneered for logarithmic theories in [42] and is now referred to as the standard 
module formalism. It has since been shown to yield correct (or, at least, sensible) results rather generally, 
see [22,43-47]. We refer to [48,49] in which general features of the standard module formalism are discussed. 

Armed with this character information, we turn to the detailed structure of the fusion product of two Kac 
modules. The additional structural data may be obtained, at least in certain examples, by explicitly constructing 
(a truncation of) the fusion product using the Nahm-Gaberdiel-Kausch fusion algorithm [28,50]. In Section 4, 
we detail an explicit example in order to illustrate how one typically employs this algorithm, noting that a priori 
knowledge of the fusion product’s character can lead to significant simplifications. We also use this example 
to emphasise that completely identifying the result, up to isomorphism, may require computing additional 
indecomposability parameters [28,51], in particular the N = 1 analogues of the logarithmic couplings of [30]. 
Moreover, we demonstrate that structural theorems for more general classes of Neveu-Schwarz modules may 
also be used to significantly simplify the identification of the fusion product. 

Section 5 then summarises the fusion rules that we have obtained for Neveu-Schwarz Kac modules by com¬ 
bining the Nahm-Gaberdiel-Kausch fusion algorithm with the Verlinde formula and Neveu-Schwarz structural 
theorems. An infinite family of fusion rules are subsequently conjectured and encoded in polynomial ring struc¬ 
tures. We conclude with a brief discussion that looks toward generalising these results to the Ramond sector and 
outlines how the formalism being developed will contribute to our long term research goals. This is followed by 
three appendices. Appendix A provides a detailed derivation of the coproduct formulae that underlie the Nahm- 
Gaberdiel-Kausch fusion algorithm and a description of the usage of the algorithm itself. While this material 
roughly follows [28,52], it makes the explicit computations in Section 4 self-contained as well as prepares the 
reader for the rather more involved generalisation, detailed in [12], needed to discuss the Ramond sector. Appen¬ 
dix B then introduces the Neveu-Schwarz staggered modules, following [48,53,54] and [14], that appear in certain 
fusion products. We prove a few basic structural results before reviewing the definition of logarithmic couplings 
and explaining how these parameters may be computed. Appendix C contains a selection of fusion products 
computed by combining the information provided by the Verlinde formula, the Nahm-Gaberdiel-Kausch fusion 
algorithm and the theory of staggered modules. 
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2. Background and conventions 


The N = 1 superconformal algebras may be defined as the Lie superalgebras spanned by the even (bosonic) 
modes L n and C, and the odd (fermionic) modes Gk , subject to the commutation relations 


[fmi -T n ] — (TO Tl)L m -\. n T 

{ Gj > Gk } = + 


4 j 2 - 1 


12 


TO — TO 

12 

^•+fc=o C. 


^m+n—()Gj [L m ,G k ] — ^ ^ TO k^jGm+k) 


( 2 . 1 ) 


[L m ,C] = [Gj,C] =0. 


In (2.1), we will take m,n £ Z and j,k stipulating that we are studying the Neveu-Schwarz algebra. 

Taking j, k £ Z results in the Ramond algebra instead. The central element C will be taken to act in all 
representations as a fixed multiple c of the identity operator, called the central charge. Formally, we thus 
consider the quotient of the universal enveloping algebra of the Neveu-Schwarz algebra by the ideal generated 
by C — c id; we will also refer to this quotient as the Neveu-Schwarz algebra. 

In field-theoretic terms, the Neveu-Schwarz superalgebra extends the Virasoro algebra by the modes of a 
fermionic primary held of conformal weight |. More precisely, the fields generated by the bosonic and fermionic 
modes are, respectively, the energy-momentum tensor T(z) and its superpartner G(z): 

T(z) = J2 L n z ~ n ~ 2 , G(z)= J2 G jZ -J~ 3 / 2 . (2.2) 

ne z je z+i/2 

The operator product expansions equivalent to (2.1) then take the form 

c/2 2 T(w) dT(w) 


T(z)G(w) 


T(z)T(w) 

[z — w) 

3 G{w)/2 dG(w) 


( z-w ) z-w 
supplemented by the locality condition T(z)G(w) = G(w)T(z). 


(z — w) 2 
G(z)G(w) 


z — w 
2c/3 


2 T(w) 
(z — wY z — w 


(2.3) 


+ 


2.1. Highest-weight modules. Neveu-Schwarz highest-weight theory works as one would expect. The trian¬ 
gular decomposition splits the superalgebra into the span of the positive modes L n and Gj, with n,j > 0, the 
negative modes L n and Gj , with n,j < 0, and the zero modes Lq and C = c id. A highest-weight vector Vh 
is therefore characterised by its conformal weight h (we regard the central charge of the module as implicitly 
fixed) and satisfies 

L n Vh = GjVh = 0 for n,j > 0; L 0 v h = hv h . (2.4) 

The Verma module Vh is then constructed from Vh as an induced module and it has a unique simple quotient 
that we shall denote by . 

These Neveu-Schwarz modules are naturally Z 2 -graded by choosing the parity of the highest-weight vector 
to be even or odd. Structurally, this choice makes no difference, but it is sometimes useful to keep the parity 
explicit. When this is the case, we affix a superscript sign ± to the module, with the sign matching the parity 
chosen for Vh- We will generalise this convention to all the indecomposable Neveu-Schwarz modules considered 
here by matching the superscript sign to the (common) parity of the vectors of minimal conformal dimension. 

The standard parametrisation suggested by the Neveu-Schwarz analogue of the Kac determinant formula is 


c = ^ -3(f+ t J ), 




- 1 


r 1 - 


rs — 1 


-t, 


(2.5) 


where t £ C\ {0}. The Verma module V r>s = Vh r s is then reducible when r and s are positive integers satisfying 
r = s mod 2. 1 If t is rational, then this parametrisation may be written in the form 


( = £ c= j 

p' 2 \ pp' 




{p'r-psf - (p' -pf 

8 pp' 


( 2 . 6 ) 


^In contrast, Verma modules for the Ramond algebra turn out to be reducible for positive integers r and s satisfying r 7 ^ s mod 2. 
However, the parametrisation for h r , s is also corrected in the Ramond sector by adding jq. 
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Point Link Chain Braid 



Figure 1. The singular vector structure, marked by black circles, of Neveu-Schwarz Verma 
modules. Arrows from one singular vector to another indicate that the latter may be obtained 
from the former by acting with a suitable polynomial in the L n and Gj. Note that t > 0 
corresponds to c < § and t < 0 corresponds to c ^ • 

where one customarily takes p = p' mod 2 and gcdjp, y (p' ~ P) } = 1- The N = 1 superconformal minimal 
models [2,34,55] are built from the Neveu-Schwarz simple liighest-weight modules £,. jS = C} lr a with 1 ^ r p— 1, 
1 ^ s ^ p' — 1 and r = s mod 2, as well as their Ramond counterparts. However, we are not studying these 
minimal models, so we do not insist, for example, that p,p' ^ 2. 

In most respects, the highest-weight theory for the Neveu-Schwarz algebra parallels that of the Virasoro 
algebra. In particular, the submodules of a Verma module are generated by singular vectors and the maximal 
dimension of the space of singular vectors of any given conformal weight is 1. The submodule structure of a 
Verma module Vh then reduces to knowing its singular vectors and here the possibilities exactly mirror those of 
the Virasoro Verma modules [7,8]. We reproduce the structures diagrannnatically in Figure 1. As mentioned 
above, if h ^ /i r , s for any positive integers r and s with r = s mod 2, then Vh is simple. On the other hand, if t 
is irrational, then the maximal proper submodule of the Vh with h = h rtS , r = s mod 2, is always simple. It is 
generated by a singular vector of depth \rs , meaning that its conformal weight is h r , s + \rs. 

We will exclusively focus on the case in which t = p/p’ is rational and positive, with h = h r ^ s and r = s mod 2. 
If r is a multiple of p, or s is a multiple of p', then the singular vector structure of Vh is represented by the 
infinite chain diagram in Figure 1. Otherwise, the structure corresponds to the infinite braid diagram; this 
latter case is the one relevant to the study of minimal models. In both cases, a singular vector is always present 
at depth ^rs, though there may be other singular vectors at other depths. 

All of this information may be conveniently summarised in the Neveu-Schwarz analogue of the extended Kac 
table, see Figure 2. This simply tabulates the values of h TjS as r and s run through the positive integers, subject 
to r = s mod 2. 2 To make contact with the above structural results, we partition the extended Kac table into 
three subsets as follows: 

• If p divides r and p 1 divides s, then we say that (r, s) is of corner type in the extended Kac table. 

• If p divides r or p' divides s, but not both, then (r, s) is said to be of boundary type. 

• If p does not divide r and p' does not divide s, then (r, s) is said to be of interior type. 

Summarising, corner and boundary type Verma modules have singular vectors arranged in chains whereas 
interior type Verma modules have a braided pattern of singular vectors. We remark that when p = 1 or p' = 1, 

2 The “gaps” in the table, when r ^ s mod 2, correspond to conformal weights of representations of the Ramond algebra. 
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there are no interior entries in the extended Kac table, and if p = p' = 1, then there will be no boundary entries 
either. We illustrate this with three pertinent examples of extended Neveu-Schwarz Kac tables in Figure 2. 


2.2. Neveu-Schwarz Fock spaces. The N = 1 superconfornral algebras have a free field realisation in terms 
of a free boson and a free fermion, the latter taken in the free fermion Neveu-Schwarz or Ranrond sector to 
obtain the corresponding superconformal sectors. In particular, the Neveu-Schwarz algebra acts on the tensor 
product of any bosonic Fock space with the vacuum fermionic Fock space. We shall refer to such tensor products 
as Neveu-Schwarz Fock spaces. 

At the level of operator product expansions, one starts with a free boson field a(z) = ’Yh ne j 1 a n z ~ n ~ 1 and a 
free fermion field b(z) = z- 1/2 satisfying 


a(z)a(w) 


1 

7 \ 2 > 

(z — w) 


b{z)b(w) 



(2.7) 


These are then used to construct the energy-momentum tensor and its superpartner: 

T{z) = i : a(z)a(z) : + Q-da(z) + ^ : db{z)b{z) : , G(z) = a(z)b(z) + Qdb(z). (2.8) 


Here, : • • • : denotes normal ordering. It is straightforward to check that these fields satisfy the operator 
product expansions (2.3) with c = | — 3 Q 2 . To match the central charge of (2.6), we set 


0 = 



(2.9) 


The Neveu-Schwarz Fock space J~\, being the tensor product of a free boson Verma module with the free 
fermion vacuum module, is generated by a highest-weight vector v\ satisfying 


a n v\ = bjv a = 0 for n,j > 0; a 0 v x = \v\. 


( 2 . 10 ) 


This generator v\ then has conformal weight 


hx = 

which will coincide with a weight h r<s 


4 pp'( x -Q / 2 ) 2 - & -p ) 2 

2 A( - ™ 8 pp' 

of the extended Kac table, given in (2.6), when 


( 2 . 11 ) 


A = \ r ,s = —a!{r — 1) + a(s — 1). 


( 2 . 12 ) 


Here, we have introduced the quantities 


a = 



a = 


and note the symmetries 


^r+p,s — ^r,s r, \!PP r i ^r,s+p' — ^ r,s T ^ \/ PP' 


'^r4-p,s+p / — ^r,s? 


(2.13) 


(2.14) 


for later use. 

Fock spaces are always simple as modules over the product of the free boson and fermion (universal enveloping) 
algebras. However, this need not remain true upon restricting to the Neveu-Schwarz algebra. Specifically, the 
Fock space J~\ will only be simple as a Neveu-Schwarz module when A ^ A rjS for any r, s € Z with r = s mod 2. 
For t € Q+, we describe the submodule structure of F\ as islands, a chain, or a braid, illustrating the possibilities 
in Figure 3. These structures mirror those of the Feigin-Fuchs modules of the Virasoro algebra. More precisely, 
when (r, s) is a corner/boundary/interior type entry of the extended Kac table, then the submodule structure 
of J> jS = J~\ r a is of islands/chain/braid type. There are two possible structures for chain and braid type 
Fock spaces F ryS . corresponding to the fact that these Neveu-Schwarz modules are not isomorphic to their 
contragredient duals Fq-\ r s = F- r - s . We shall describe the detailed structure of the F r s in the rest of this 
section. 
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Figure 2. Parts of three of the extended Neveu-Schwarz Kac tables for c = |, c = — 
c = 0. The rows of the tables are labelled by r = 1, 2,3,... and the columns by s = 1,2, 
Interior points are shaded dark grey, boundary points are shaded light grey, while corner points 
are white. 
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Point Islands 


Chain 


Braid 



Figure 3. The structure of Fock spaces over the Neveu-Schwarz algebra, for t £ Q + (c ^ |). 
Each black circle represents a subsingular vector and the arrows represent the action of the 
algebra as in Figure 1. The two braid diagrams are actually reflections of one another, but the 
repetition reminds us that braid type Fock spaces are not self-contragredient. 


It is important to note that, unlike the Verma modules of the previous section, the submodules of the Neveu- 
Schwarz Fock spaces are not necessarily generated by singular vectors. Instead, submodules may be associated 
with subsingular vectors, these being vectors which become singular in an appropriate quotient of the parent 
module. Unlike singular vectors, the Fock space subsingular vectors are not usually unique up to normalisation 
because subsingularity is preserved by adding any element of the submodule by which one quotients to obtain a 
singular vector. In Figure 3, we may regard each black circle as representing a subsingular vector, modulo the 
aforementioned non-uniqueness, generating submodules of the Fock space. 

It is natural to associate subsingular vectors with the simple quotient of the submodule that they generate. 
The circles in Figure 3 thus also represent the (simple) composition factors of the Fock space. For point and 
islands type Fock spaces, there are no arrows in the associated diagrams, indicating that these Fock spaces are 
semisimple. For chain and braid type Fock spaces, the composition factors may be partitioned into three classes 
according as to whether the arrows incident on the corresponding circle are all pointing towards it, all pointing 
away from it, or there are some pointing towards and some pointing away. The latter class is empty for chain- 
type modules. When all arrows point towards the circle, it represents a singular vector generating a simple 
submodule. The direct sum of the composition factors in this class therefore gives the maximal senrisinrple 
submodule, also known as the socle of the Fock space. Similarly, the maximal semisimple quotient, also called 
the head of the Fock space, is the direct sum of the composition factors in the class corresponding to all arrows 
pointing away. 

The conformal weights of the subsingular vectors of the Fock spaces J> iS , with r,s£ Z + and r = s mod 2, 
coincide with those of the singular vectors in the Neveu-Schwarz module V r . s . It therefore follows that T r ^ s has 
a subsingular vector of depth ^rs and that its contragredient dual T- r - s does too. One can verify that this 
subsingular vector is always associated to either the socle or the head of the Fock space (its circle in Figure 3 has 
either all arrows pointing towards it or all arrows pointing away from it, respectively). With the parametrisation 
(2.12) that we chose above, it turns out that the depth \rs subsingular vector is always associated with the 
head for r, s £ Z + and with the socle for r, s £ Z_. This structural identification may be extended to all r, s £ Z 
by using the symmetries (2.14). 

This realisation fixes the structure of a chain type Fock space J> jS uniquely. One only has to determine 
whether the depth ^rs subsingular vector belongs to the socle or the head and find the number of subsingular 
vectors of lesser depth; this is sufficient to distinguish between the two possibilities in Figure 3. When J> iS is 
of braid type, this information is not quite sufficient. At every other horizontal level in the braid type pictures 
in Figure 3, there is one singular and one (non-singular) subsingular vector (at the other horizontal levels, 
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excepting the highest, no vector is singular). To identify which is which, given their depths, the following fact 
is germane: if the depth of the singular vector at a given horizontal level is greater than that of the subsingular 
vector at the same level, then it will also be greater at the other horizontal levels (and vice versa). One may then 
check which has greater depth in a given module because one knows the nature of the depth ^rs subsingular 
vector. 

2.3. Kac modules. The reducible Fock spaces J> >s , with r, s £ Z and r = s mod 2, are not themselves the 
modules of central interest here. Rather, it is a certain related class of modules that we shall call (Neveu- 
Schwarz) Kac modules, following [20-22], that take centre stage. These are indexed by positive integers r and s. 
Over the Virasoro algebra, Kac modules were introduced non-constructively in [23,26,27,56] in order to describe 
the boundary sectors of the scaling limits of certain integrable lattice models. Their characters were determined 
in many examples and the results were compatible with Kac modules being identified with certain quotients of 
Verrna modules. However, more recent consistency checks have led to a different proposal [20] for the identity 
of Virasoro Kac modules, at least for some models, as submodules of Fock spaces rather than as quotients of 
Verrna modules. This proposal has been generalised in [22] where evidence verifying this identification in many 
non-trivial examples is presented. 

Neveu-Schwarz Kac modules were recently considered from a lattice point of view in [14]. Although this 
analysis only studied the action of To on certain examples and so obtained only a bare minimum of structural 
information, it is reasonable to expect that these Neveu-Schwarz Kac modules may likewise be identified with 
submodules of Neveu-Schwarz Fock spaces. Extrapolating the results of [22] leads us to define the Neveu- 
Schwarz Kac module /C r , s , with r,s € Z + and r = s mod 2, to be the submodule of the Neveu-Schwarz Fock 
space J> iS that is generated by the subsingular vectors of depth strictly less than |rs. Note that this does not 
exclude the possibility that /C r , s may have a singular vector of depth greater than | rs. 

We illustrate this definition with examples of Kac modules for (p,p') = (2,4) (c = 0). First, we determine 
the structure of the corresponding Fock spaces in Figure 4. Remembering that J> iS = J-' r + Py s+p', by the 
symmetries (2.14), we arrive at the Kac module structures depicted in Figure 5. For general p and p ', with 
gcdjp, Tj(p' — p)} = 1, the structure of the Kac module /C rtS is indicated, for small r and s, in Figure 6. 

3. A Verlinde formula 

For rational theories, a relatively efficient route to computing the fusion rules is to first determine the 
modular S-transforms of the characters of the simple modules and then apply the Verlinde formula. For 
logarithmic theories, such a formula cannot compute the fusion multiplicities themselves, because characters 
cannot distinguish between an indecomposable module and the direct sum of its composition factors. Instead, 
it is natural to require that a logarithmic Verlinde formula computes the structure constants of the projection 
of the fusion ring onto an appropriate ring of characters. 3 In many cases, this character ring may be identified 
with the Grothendieck ring of fusion, in which indecomposable modules are identified with the sum of their 
composition factors. 4 The expectation is therefore that the Verlinde formula determines the character, or 
equivalently the composition factors, of a fusion product. 

There are a few provisos to this claim. Technically, a Grothendieck ring of fusion will only exist if fusing with 
any given Neveu-Schwarz module defines an exact functor from our chosen category of Neveu-Schwarz modules 
to itself. While this is not true for general Neveu-Schwarz modules (the first example over the Virasoro algebra 
was given in [33]), the modules for which it is true form a subring of the fusion ring [64] (assuming that fusion 
defines a tensor structure on the appropriate category of Neveu-Schwarz modules). It has been conjectured that 
modules defining boundary sectors in a consistent boundary conformal field theory do define exact functors 
under fusing. We will therefore assume that the Kac modules have this property; consequently, this property 

3 Other, typically model-dependent, Verlinde-like formulae have been proposed for logarithmic conformal field theories in [57—63]. 
3 This identification assumes that the characters of the simple modules are linearly independent. This is indeed the case here, 
provided that one does not distinguish modules if they only differ through the parity of their highest-weight vector. If one does 
want to distinguish these modules, then supercharacters should be used as well as characters. In the latter setup, the sign with 
which a character appears becomes relevant. 
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Figure 5. Examples of the structures of Neveu-Schwarz Kac modules when p = 2 and p' = 4. 
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Figure 6. A depiction of the structures of the Kac modules /C, %s as (r, s) varies over (a part of) 
the extended Kac table. The genuine Kac table, bounded byl<r<p—1 and 1 < s < p' — 1, 
is represented by the dark grey rectangle in the upper-left corner. These are interior points of 
the extended Kac table and the light grey and white areas correspond to boundary and corner 
type labels as in Figure 2. If p = 1 or p' = 1 (or both), then the possible structures correspond 
to removing the rows or columns (or both) that contain interior labels. 


will be shared by the modules that are generated from the Kac modules by fusion. Thus, the subring of the 
Neveu-Schwarz fusion ring that is generated by the Kac modules is hereby assumed to possess a well defined 
Grothendieck ring of fusion. This subring is central to our investigations and we shall begin our analysis by 
exploring the detailed structure of its Grothendieck ring. 


3.1. Modular transformations. We start with the study of the modular transformation properties of Neveu- 
Schwarz characters. Those of the Fock spaces are particularly accessible: 


ch[j- A ](r)=g^- c / 24 n 

i=i 


1 + gj-l/2 _ g (A-Q/2) 2 /2 /tf 3 (l -q) 


1 — qi 


v(q) 


v(q) 


(q = e 27r,T ) 


(3.1) 


Here, h\ denotes the minimal conformal weight of J~\ and we employ the Dedekind eta and Jacobi theta 
functions 

oo oo 

v(q) = q 1/24 II i 1 ~ <?) ’ ^ (z; q) = n i 1 + zq j ~ i/2 ) (i - q 3 ) (i + V" 1/2 ) • ( 3 . 2 ) 

i=i 

Note that the character formula (3.1) does not distinguish between F\ and its contragredient dual Jq_a- 
The modular S-transformation of these characters is then 

pOO 

ch[J- A ] (—1/r) = S[F X ^ J- M ]ch[j- M ] (r) d/x, S[T X —> = 2cos[2tt(A - Q/2)(/x - Q/2)], (3.3) 

Jq/2 


where we restrict the range of the integral to [Q/2, oo) in order to integrate over linearly independent characters. 
This may be easily verified using a standard gaussian integral, convergent for Imr > 0. However, because this 
S-transform requires a continuous family of Fock space characters, we expect to encounter singular distributions 
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at some point. As it can be confusing to allow endpoints to the integration domain when computing with these 
generalised functions, we will redefine the above S-transformation, once and for all, so that the integration range 
is open: 

/ OO 

S [Fx -A ch [^] (r) dp, S [Ja -a F m ] = cos [2 tt(A - Q/2){p- Q/2 )\. (3.4) 

-oo 

The price to pay for extending to all /i £ R is that we must now identify T p and its contragredient in all 

computations involving characters, in particular when we employ the Verlinde formula. 

We remark that the kernel S \iF\ -A of this S-transform, the analogue of the S-matrix for rational theories, 
is symmetric (S \T\ -A J 7 ^] = S [j 7 ^ -A _Fa] ) and unitary: 


/ oo />oo 

S[F X -A -A J^]* dp = / cos[27 t(A - Q/2)(p - Q/2)] cos[2tt(u - Q/2)(/z - Q/2)] d// 

-oo J — OO 

1 /*°° 

= - / (cos[27r(A — u)/t] + cos[27t(A + u — Q)/i]) d^t 

^ J —oo 


= - \8{v = A) + 8{v = Q - A)] = 8{v = A). 


(3.5) 


Here, the final equality is justified by the fact that we must identify A and Q — A in such computations. Because 
the S-transformation kernel is real, this computation also proves that S squares to the identity operator. In 
other words, conjugation is trivial in this theory (at the level of characters), as one would expect. 

These results accord with our general expectations for good modular properties and indicate that we may 
expect meaningful results from the Verlinde formula. Indeed, in the general formalism for the modularity of 
logarithmic conformal field theories proposed in [48] and refined in [49], the T\ may be regarded as the standard 
modules : They are simple for almost all A € R, with respect to the Lebesgue measure used in (3.4); the simple 
F\ are termed the typical modules. The atypical modules are then those corresponding to A = A rjS , for r,s € Z. 
In particular, the J~ r s are not simple, hence are atypical. 

The Kac modules K. r>s of Section 2.3 are therefore also atypical and their characters are given by 


ch[/C r , s ] = chfj^J - ch[j 7 A_ r s . 


(3.6) 


We immediately obtain their S-transformations: 

/ OO 

S[/C r ,g -A- ch [J 7 ^] (r) dp, 

-oo (3.7) 

S[/C„g -A Ffj] =S[J 7 A r s Sf.7 7 \_ r s a/J = 2sin[27rra'(/z - Q/2)] sm[2nsa(p - Q/2)]. 

Note that the result of S-transforming the character of a Kac module is an integral over the Fock space characters. 
As the Fock spaces are the standard modules of the theory, their characters give the canonical topological basis 
in which to express all characters. 5 In particular, it is not clear that the quantity S[/C r . s —> K, r ',s’\ is well 
defined. We will therefore perform all subsequent computations in the basis of standard characters without 
further comment. 

Finally, we recall that Kac modules were only defined for r,s £ Z+. If the character formula (3.6) is extended 
to general r, s € Z, then we obtain 

ch[/C_ riS ] = — ch[/C r>s ] = ch[/C r) _g], ch[/C r> o] = ch[/Co,g] =0, ch[/C_ rj _ s ] =ch[/C r>s ]. (3.8) 


These formulae are important for interpreting the results of general Verlinde computations. 


3.2. Verlinde products. Assuming, as discussed above, that the fusion product “x” descends to a well defined 
product ‘W on the Grothendieck ring of characters, we may decompose a character product into a linear 


®We refer to this basis as topological because one might only recover a given character, for example that of a simple atypical module 
C r ,s, as an infinite sum of Fock space characters. In this case, the convergence is that of formal power series — the contribution to 
the multiplicity of a given weight space is zero for all but finitely many terms in the sum. 
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combination of Fock space characters: 

ch[A/l x Af] = ch[Ad] El ch[A f] = 


dv - 


(3.9) 


Here, Ad and A f are Neveu-Schwarz modules and we recall that the Fock space characters are the preferred 
(topological) basis for the space of all Neveu-Schwarz characters. The multiplicities Nare the Verlinde 
coefficients and are computed, in terms of the S-transformation kernel, by the Verlinde formula: 




S [M -a T p \ S [AT -a T p \ S [T v -a T p ] ’ 

s \JC i,i —a Tp 


dp. 


(3.10) 


Here, /Cyi plays the role of the vacuum module. Note that this Kac module is generated by a higliest-weight 
vector of conformal weight 0 that is annihilated by both L_i and G_ 1 / 2 . 

It is now straight-forward to compute the character of fusion products. The unitarity of the S-transformation 
implies that the vacuum module K .yi is the unit of the character product: 


ch[/Ci,i x A f\= ch[/Ci,i] H ch[A f\ = ch[AA]. 

A somewhat less trivial example involves the fusion of /C 34 with an arbitrary Fock space: 


(3.11) 




s [/C 3i i 


— OO 
nOO 


s [/Cl,l 


T b 


dp 


sin[67ra'p] 


cos[27r(^ — Q/2)p\ cos[27t(u — Q/2)p] dp 


sin[27r a'p\ 

OO 

(1 + 2cos[47ra'p])(cos[27r(/x — v)p] + cos[27 i{p + v — Q)p ]) dp 


= 8{v = p — 2a) + = p) + 5{v = p + 2a) 


=> cli[/C 3 > i] Klchfjyj 
Similar computations result in 


= Cll[j) t _ 2 a'] + Cllfjyj] + chfj) 


fi+2a' \ 


ch[/Ci >3 ] lElchfj^] =ch[J M _ 2a] + ch + chJ> M+2a ], 

ch[/C 2 , 2 ] ^ ch. \^J ~— ch [T -J- ch \p~ u,—a'+a] -\- cha] -|-ch|^^* w+ck'+q;] • 
Using Equation (3.6), we obtain the corresponding products with arbitrary Kac modules: 
ch[/C 3)1 ] Kch[/C rjS ] = ch[/C r _ 2 , s ] + ch[/C r)S ] + ch[/C r+2 , s ], 
ch[/Ci )3 ] Kch[/C rjS ] = ch[/C rjS _ 2 ] +ch[/C r , s ] + ch[/C r>s+2 ], 

ch[/C 2)2 ] Kch[/C rjS ] = ch[/C r _i >s _i] + ch[/C r _i, a+ i] + ch[/C r+l,s—l] + ch[AC r +i )a +i] . 


(3.12a) 


(3.12b) 

(3.12c) 

(3.13a) 

(3.13b) 

(3.13c) 


Here, we must employ (3.8) if the labels on the Kac modules of the right-hand side are not positive integers. 

It follows from these character products that the Kac characters span a unital subring of the Grothendieck 
ring of Neveu-Schwarz characters and that this subring is generated by ch[/C 3 > i], ch[/C 2 , 2 ] and ch[/Ci >3 ]. 6 As¬ 
sociativity then leads to an explicit general formula for Kac character products: 

r-\-r'— 1 s+s / —1 

ch[/C r>s ] E3ch[/C r q a /] = £' £' ch [JC r »,s»]- (3-14) 

r"—\i — r' 1+1 s // = |s—s'|+l 

Here, the primed sums indicate that the summation variable increases in steps of two. We mention, for later 
purposes, the following special case: 


ch[/C r> i Elch[/Ci jS = ch[/C rjS . 


(3.15) 


^There is a simple exception when c = (p = p' = 1 ) because then ch[/C3 5 i] = ch[/Ci,3] and ch[/C2,2] = ch[/Ci } i] + ch[/Ci,3]. 
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4. An explicit fusion product 


In this section, we use an example to illustrate the steps involved in completely decomposing a fusion product 
and identifying its (indecomposable) direct summands. To construct the fusion product itself, we utilise the 
Nahm-Gaberdiel-Kauscli fusion algorithm, referring to Appendix A for further details concerning this technology. 
The example is fairly involved and we have chosen it in order to illustrate a wide variety of the features and 
methods that we employ to analyse more general fusion rules. Some rather more simple arguments are presented 
in Section 5.1 (though stripped of the explicit Nahm-Gaberdiel-Kausch computations). 

As our example, we consider the fusion of the Neveu-Schwarz Kac module K,-\ .3 with itself at central charge 
c = 0 (p = 2 and p' = 4). A part of the extended Neveu-Schwarz Kac table for c = 0 appears in Figure 2. We 
remark that /C 1 , 3 , unlike most Kac modules, is a highest-weight module; indeed, it is generated by a highest- 
weight vector v of conformal weight = 0. We may therefore identify AC 1,3 as the quotient of the Verma 
module Vo by the submodule generated by the singular vector of conformal weight |. Thus, 


(L-iG-1/2 - 



v = 0 


(4.1) 


in AC 1 , 3 . For simplicity, we shall assume throughout that v is even. 

First, we determine the character of the fusion product using the Verlinde formula. Specifically, Equa¬ 
tion (3.13b) gives 

ch[/C 1)3 ] Kl ch[/C 1>3 ] = ch[/C M ] + ch[/Ci )3 ] + ch[/Ci, 5 ]. (4.2) 

However, each of the Kac modules appearing on the right-hand side is reducible, with two (simple) composition 
factors each, so we learn that the fusion product has six composition factors in all: 


ch[/Ci, 3 x /Ci,3 — 2 ch[£o] + 2 ch1 /2] T ch[>^3/2] + ch^.5]. 


(4.3) 


Here, Ch denotes the simple highest-weight module whose highest-weight vector has conformal weight h, as in 
Section 2.1. To understand how these six simple modules are glued together to form the fusion product, we will 
partially construct the product module and explicitly analyse the action of the Neveu-Schwarz algebra upon it. 

To construct the fusion product of AC 1,3 with itself, we first calculate its special subspace. This is defined (see 
Appendix A.2) to be the (vector space) quotient of ACi ,3 by the action of the algebra generated by the Virasoro 
modes L n , with n ^ —2, and superfield modes Gj, with j ^ — |, leaving only linear combinations of vectors in 
which L_ 1 and G_i / 2 act on v. Imposing the singular vector relation (4.1), we find that L_\G _\i 2 v = \G_^/ 2 v 
must be set to 0 in the special subspace; generalising this shows that the special subspace is three-dimensional: 

ICf 3 = span{ v, G_ 1/2 v, L_iu}. (4.4) 


We will first determine the depth 0 truncation of the fusion product ACi ,3 x ACi, 3 . Naturally enough, this 
requires the depth 0 truncation of ACi, 3 - This subspace is obtained by quotienting by the action of all Neveu- 
Schwarz monomials with negative indices: 

AC° 3 = span{u}. (4.5) 

The depth 0 truncated fusion product is constructed within the tensor product of these two spaces. Thus, 

[/Ci ,3 x /Ci, 3 ]° C /Ci ® 3 0c /C i , 3 = span{u 0 v, L_iv 0 v | G _\/ 2 v 0 v}. (4.6) 

The basis vectors here have been partitioned into even and odd parities, using a vertical delimiter, recalling that 
v has been assumed to be even. To determine which subspace of this three-dimensional tensor product is the 
depth 0 fusion product, we search for spurious states. These are linear dependence relations (see Appendix A.2) 
that may be derived in ACg s 3 0c /Cj 3 when we impose the Neveu-Schwarz algebra action defined by the fusion 
coproduct formulae. Inspection of the composition factors (4.3) shows that there must be two vectors of 
conformal weight 0 in the depth 0 product, hence there can be at most one spurious state. 
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We search for spurious states by implementing the singular vector relation (4.1). This will require the 


following cases of the master formulae (A. 18) derived in Appendix A.l: 

A(G_]y 2 ) = G—i /2 0 1 + Mi 1 ® G_i/ 2 i (4.7a) 

A(L_i) =L_i<g>l + 10 L_i, (4.7b) 

A(G_ 3 / 2 ) = G_ 1/2 0 1 — • • • + Mi 1 0 G_ 3 / 2 , (4.7c) 

G_3/2 0 1 = A(G_3/2) + • • • + mi [l 0 G_ 1/2 + ■■•]. (4.7d) 


Here, A is the fusion coproduct, m 1 = ±1 is the parity of w\ when the formula is applied to w\ 0 W 2 , and the 
dots stand for infinite numbers of omitted terms which will not contribute to this calculation. First, we note 
that all Virasoro and Neveu-Schwarz modes, except Lq, will act as the zero operator on a depth 0 space. In 
particular, A(G_i/ 2 ) = A(£_i) = A(G_ 3 / 2 ) = 0 , so that 

0 = A(G_ 3 / 2 )u 0 V = G_ 1 / 2 V 0B + D0 G_ 3 / 2 f = G_ 1 / 2 V 0t) + 2»0 L- 1 G-I/ 2 V 

— G_i/2V 0 V 2 L — \V 0 G_iJ2V — G_ij2V 0D + 2 G _\j2^ J — 1+ 0 V = G _\j2V 0^42 L—\G_\/2V 0 V 

— G_ij2V 0 V + G_3/2^ 0 U = G_iJ2V 0 V 4 U 0 G_i/2V — 0. (4.8) 

In this calculation, we have used (4.7c), (4.1), (4.7b), (4.7a), then the commutation relations (2.1), (4.1) again, 
(4.7d), and finally (4.7a) again. We have also assumed that the highest-weight vector v has even parity. In any 
case, the right-hand side of (4.8) is identically zero which means that we have failed to find a spurious state. 
Replacing v 0 v by G_ 1 / 2 V 0 v or L_\v 0 v in this calculation likewise fails to uncover any spurious states. 

We therefore assert that there are no spurious states to find and that the depth 0 fusion product is three- 
dimensional. It only remains to determine the action of Lq, that of the other modes being trivial. To this end, 
we need three additional auxiliary formulae: 


A(l/ 0 ) — L —1 0 1 + Lq 0 1 + 10 Lq, (4.9a) 

L— 2 0 1 = A(L_ 2 ) + • • • + 1 0 L —1 10 Lq + • • • , (4.9b) 

LtiV = --G_ 3 / 2 G_i/ 2 u + L- 2 v. (4.9c) 

The last is a consequence of the singular vector relation (4.1). The action of Lq is now given by 


A(Lo)i ) 0 v = L^iv 0 v, (4.10a) 

A(Lo)G_i/ 2 u 0 » = L_iG_i/ 2 u 0 » + \g_ i/ 2 v ®v= iG_ 3 / 2 u 0 t) + ^G_ i/ 2 u 0 v 

= 0 G_i/ 2 r + ^G_i/ 2 r 0 u = 0, (4.10b) 

A(Lo) L_\v 0 b = L 2 _ x v 0 » + L-\v 0 » = —-G_ 3 / 2 G_i/ 2 f 0 » + T_ 2 V 0 v + L _\v 0 

= —G_ 1 / 2 V 0 G— 1 / 2 V + t)0 L_iv + L_iu 0v = —G'^_ 1 / 2 'C 0v = —L-iv 0 u. (4.10c) 

With respect to the ordered basis (4.6), we have 


A (L 0 ) = 


/ 0 0 

°\ 

1 h 

A 

\o 0 

0 / 


(4.11) 


where we have partitioned the matrix to indicate the separation into even and odd basis elements. We conclude 
that the depth 0 fusion product is spanned by two vectors of conformal weight 0 , one even and one odd, and 
one even vector of weight 4. 

The depth 0 result therefore accounts for three of the six composition factors of the fusion product, namely 
both of the Cq factors and one of the £ 1/2 factors. The remaining factors, £ 1 / 2 , £ 3/2 an d £ 5 , must appear as 
descendants of these via the action of the negative modes; otherwise, they would have appeared in the depth 0 
calculation. The factor £ 1/2 can only descend from one of the £q factors, but once this is fixed there are still 
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three consistent possibilities, ignoring parities, for identifying £ 3 / 2 and £5 as descendants: 


0 : 


1. 
2 • 


3. 

2‘ 






5: 


















(4.12) 


To distinguish between them, we must construct the fusion product to greater depth. 

We therefore turn to the depth 4 calculation in which vectors are set to zero if they may be obtained from 
other vectors by acting with linear combinations of Neveu-Sclrwarz monomials whose indices are negative and 
sum to at most —1. The special subspace of K. does not change, but now we consider its depth 4 truncation 
which is spanned by v and Thus, the depth 4 fusion product will be contained within a six-dimensional 

space: 


[£ 1,3 x C span{u ® v, \v ® v, G_i/ 2 v ® G_ \i 2 v I G_i/ 2 v (g> v, v <8> G_x/ 2 v, L_iv ® G_i/ 2 v }. (4.13) 


Comparing with the three possible structures (4.12) for the fusion product, we see that allowing descendants 
by G_ 1 / 2 , as well as the depth zero vectors, always leads to five depth 4 vectors with conformal weights 0, 0, 
\ and 1. This indicates that there is precisely one spurious state to find. 

The calculation proceeds in much the same manner as before. The difference is that because we are computing 
to depth 2 , we may no longer assert that A(G_i/ 2 ) = 0 (nor that A(G + i/ 2 ) = 0). Using (4.7c), (4.1) and 
(4.7b), we quickly arrive at 


0 — A(G_3/2) , e v — G_.i j 2 v ® v T v C) G_ 3/ 2 v — Gj 2 v ® v T 2 v L—\G_\j 2 v 
= G^x/ 2 v (8 v — 2 L_\v ® G_i/ 2 v. (4-14) 


The right-hand side has been expressed in terms of the basis elements (4.13) and the fact that it does not 
vanish identically means that we have found a spurious state. More precisely, it means that this relation must 
be imposed in the depth 4 fusion product. We have searched for more independent spurious states, but found 
none in accord with the structural arguments above. 

Imposing this relation reduces the dimension of the depth 4 fusion product from 6 to 5. Computing the 
action of Lq on this space is now straight-forward. With respect to the ordered basis consisting of the first five 
elements of the right-hand side of (4.13), we obtain 


A(Lq) = 


/o 0 0 

0 

0 


/o 0 0 

0 o\ 

1 0 -4 

0 0 


1—1 

rH|(M 

0 

0 0 

0 4 1 

0 0 


1-Hl 

0 

0 

0 0 

0 0 0 

1 1 

2 2 


0 0 0 

0 0 

0 

0 

^0^ 

1 i / 

2 2/ 


\0 0 0 

0 1) 


(4.15) 


where we also indicate the Jordan canonical form. While we do find the expected conformal weights, more 
interesting is the presence of a rank 2 Jordan block for the weight indicating the presence of a staggered 
submodule in the fusion product (see Appendix B). 

To determine which of the three possibilities of (4.12) is realised by the fusion product, we can repeat the 
above computations to depth 4 and show that the fusion product has no submodule isomorphic to Cq, that is 
that no weight 0 vector is annihilated by both singular combinations G_ 1/2 and L_iG_i /2 — ^G_ 3 / 2 - This fact 
implies that the fusion product corresponds to the leftmost possibility in (4.12). We mention that the required 
computation is rather tedious by hand, involving one spurious state in a 12 -dinrensional space, but is practically 
instantaneous in our computer algebra implementation. 
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The full structure of the fusion product is 



Am, 3 x A, 3 — /Ci,i ® S 1 


0,1 

4) 


(4.16) 


where S® 4 denotes a staggered module described by the short exact sequence 


/C 


1,3 


c o,i 

D l,4 


/C 


1,5 


0 . 


(4.17) 


The notation here derives from (4.17) in that £44 has a submodule isomorphic to the Kac module with labels 
(r, s) = (1,4) — (0,1) = (1,3) and the quotient by this submodule is isomorphic to the Kac module with labels 
(r, s) = (1,4) + (0,1) = (1,5). We will use the obvious extension of this notation to describe more general 
staggered modules in what follows (see Appendix B.l ). 7 

With our depth f computation, we are now able to check every aspect of (4.16) except for explicitly verifying 
the arrow from the subsingular vector of conformal weight 5 to the Lo-eigenvector of weight 4. Unfortunately, 
this would require computing fusion truncations to depth 5 which is well beyond the current limits of our 
computer. Below, we will discuss an alternative means of checking that this arrow is present. 

However, having determined that the fusion product involves a staggered module <Sq 4 , we have to determine 
if the structure depicted in (4.16), equivalently if the short exact sequence (4.17), completely specifies its 
isomorphism class. The general theory states that this isomorphism class is characterised by its logarithmic 
coupling /?4 4 which may be determined, in this example, as follows (see Appendix B.2 for generalities): Let x £ 
544 have conformal weight 0, so that G_ 1 / 2 X is singular. Choose any y £ £4 4 satisfying (L 0 — \)y = G_ 1 / 2 X 
and note that G\j 2 y must be proportional to x. The constant of proportionality is /3° 4. 

We may compute /J 44 within the depth 4 fusion product by computing A(G_ 4 / 2 ) and A(G' + i/ 2 )- In the 
basis consisting of the first five elements of the right-hand side of (4.13), we find that 


A(G_4/ 2 ) 


/o 0 0 

0 0 \ 


/o 0 0 

o o\ 

0 0 0 

1 -1 


0 0 0 

1 0 

0 0 0 

-1 1 

, A(G +1 / 2 ) — 

0 0 0 

0 1 

1 1 1 

0 0 


1 1 1 

0 0 

Vo i i 

0 0 ) 


Vo i 0 

0 0/ 


(4.18) 


The element x £ £44 may be identified (in the depth 4 product) with the vector (0, 0, 0 11, — 1) T . We then 
solve ( L 0 — i)y = G_i/ 2 x , giving y = (0, —2, —2 | 0,0) T modulo arbitrary multiples of G_ i/ 2 x and so Gi/ 2 y = 
(0,0,0 | — 3, —1) T = —x — 2z, where z = (0,0, 0 11,1) T has conformal weight 1. Now, this appears to contradict 
the fact that Gi/ 2 y must be proportional to x. This is down to a subtlety with the computation of A(G 4 / 2 )- 
This mode should not be regarded as mapping the depth 4 fusion product into itself, but rather as a map from 
the depth ^ product into the depth 0 product. The vector z, being of conformal weight 1 and thus not in the 
depth 0 product, should therefore be set to 0 in order to arrive at the correct result: Gi/ 2 y = —x. We therefore 
conclude that (3 ®= — 1 and identify the fusion product as 


/Ci )3 x/Ci , 3 =/Ci,i®<Si°; 4 1 (-l). 


(4.19) 


Actually, we can refine this even further by keeping track of parities. We assumed in our computations that 
the minimal conformal weight vectors of both copies of /Ci ? 3 were even. The same is true for the summand /Ci ? i 
found above, though the vector x of minimal conformal weight in <Si 4 (—1) is odd. A maximally precise version 


^We emphasise that this notation differs from a similar notation 7 Zr’s that has been used to indicate certain modules over the 
Virasoro [26, 27] and Neveu-Schwarz [14] algebras. These modules are believed to arise in the continuum scaling limit of certain 
statistical models via a lattice fusion prescription and are conjectured to have Jordan blocks for Lq of rank 2, if exactly one of a 
and b is non-zero, and rank 3, if both a and b are non-zero. If the rank-2 module TZr's is staggered, then the two notations are 
believed to agree: "R-r’s — Sr’s • 
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of (4.17) and (4.19) is therefore 


/C+x/C+=/C+©<S 1 0 ; 4 1 (-l)- 


KL 


0-1 r 


*2. 


0. 


(4.20) 


We have also confirmed the logarithmic coupling /3 14 = —1 using the method described in [31]. This 
succeeds because, in this case, the exact sequence (4.17) fixes the isomorphism class of the staggered module 
<S 4 4 completely. We omit the rather tedious calculations. Instead, we indicate how to confirm this value using 
the heuristic formula (B.5) developed in [65]. For this, we first perturb the parameter t, hence the central charge 
c and Kac weights h r>s , as in (B.4) to t(e) =t + e= \ + e. We then compute the scalar product 


(x(£),G 1/2 G_i /2 x(£)) = (x(£),2 L 0 x(e)) = 2 h lt3 (e) = 2e 


(4.21) 


in the (poorly characterised) perturbed theory in which x(e) is a highest-weight vector of conformal weight 
h lt3 (e) = £. Substituting into (B.5), we obtain 




0,1 

1,4 


8t 2 2 e _ 

0 - (5 2 - 3 2 )f 2 e™ T “ 


(4.22) 


in agreement with the explicit depth 4 fusion construction. 

Finally, we point out that one can also arrive at the leftmost possibility in (4.12), without performing depth 
| calculations, by instead appealing to the theory of staggered modules (Appendix B). This is generally far 
more efficient than explicitly constructing the truncated fusion product, especially when the depth required for 
a complete identification becomes large. First, there are two independent vectors of conformal weight 4 and we 
know that only one, w say, is descended from a weight 0 vector. As the weight 0 vectors are eigenvectors of To, 
so is w. Thus, w is the Lo-eigenvector in the Jordan block at weight Now, choose a Jordan partner y for 
w so that (To — — w. Then, Proposition B.3 shows that w cannot have a singular descendant of weight 5 

unless this descendant also has a Jordan partner (take U there to be the singular combination of weight | that 
annihilates iry). This rules out the rightmost diagram in (4.12) — the corresponding staggered module simply 
does not exist. 

One can similarly rule out the middle diagram using some deeper structural results for staggered modules. 
In this case, the purported staggered module S(— 1) would be described by the following exact sequence: 


Vo 


V 3 + V 5 


S(-l) 


Vl/2 

v 3 


0 . 


(4.23) 


Here, we indicate the required logarithmic coupling which was obtained from a depth 4 computation. If we 
replace the third module in this sequence by its Vernra cover V 3 / 2 , then a staggered module S(— 1)' with this 
new exact sequence may be shown to exist (and be unique up to isomorphism) using the same methods that 
were employed in [54] for Virasoro staggered modules. Moreover, a staggered module <S(— 1) with sequence 
(4.23) will exist if and only if there exists a singular vector of conformal weight 3 in <S(— 1)'. It is not difficult 
to check this explicitly with a computer implementation of <S(—1/ — the result is that no such singular vector 
exists, hence that <S(—1) does not exist either. This rules out the middle diagram. 

A similar calculation may be used to verify the presence of the arrow in (4.16) from the conformal weight 
5 subsingular vector to the weight 4 singular vector G_i/ 2 x. If this arrow were not present, then this weight 
5 vector would have to be singular in <S 4 4 (—1). Again, an explicit search for a singular vector of this weight 
leads to no solutions, thus verifying the arrow. In principle, the arrow could instead point to the weight 0 
singular vector x £ <S 4 4 (— 1). However, this is easy to rule out because all of the positive weight vectors in the 
submodule generated by x actually belong to that generated by G_i/ 2 x. More generally, we may appeal to 
the Neveu-Schwarz generalisation of the Projection Lemma [54, Lem. 5.1] to identify the targets of such arrows 
(assuming we have shown that said arrows exist). 

To summarise, we have completely identified the fusion product considered above by combining informa¬ 
tion obtained from four distinct sources. First, the Verlinde formula (3.10) decomposed the corresponding 
Grothendieck fusion product, giving us the composition factors of the fusion product. Second, applying the 
Nahm-Gaberdiel-Kauscli algorithm to depth 0 indicated which composition factors were descendants of others 
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and which were not, while the depth 4 computation uncovered a rank 2 Jordan block for the action of Lq and 
determined the corresponding logarithmic coupling. We also determined the parities of the indecomposable 
direct summands of the fusion product from these computations. Third, the structure theory for staggered 
modules allowed us to completely fix the rest of the structure of the fusion product, except for one arrow in 
our diagrammatic representation of the structure that may or may not have been present. Fourth, the presence 
of this arrow was confirmed by explicitly showing the non-existence of a singular vector, of the appropriate 
conformal weight, in the staggered module. 

The logic that led us to the structure of this fusion product is fairly typical. We have employed it to analyse 
many further examples of Neveu-Schwarz fusion products. The conclusions that we have drawn from these 
analyses are reported in the next section. 


5. Results 

In this section, we summarise the results that we have obtained by combining the character product rules 
(3.14) with explicit Nahm-Gaberdiel-Kausch fusion computations and the structure theory of staggered modules. 
As was explained in the previous section, this combination allows us to significantly reduce the depth to which 
the fusion algorithm must be applied in order to completely identify the product. For brevity, we have only 
considered fusion rules between Kac modules, restricting to the central charges c = |, — —fjj, 0, —^ and 

corresponding to (p,p') = (1,1), (1,3), (1,5), (2,4), (2,8) and (3,5), respectively. The results obtained suggest 
conjectures for certain classes of general Kac fusion rules which we describe below. 

5.1. Fusing /C r> i with /Ci, s . Perhaps the simplest Kac module fusion products are those involving a “first 
row” module and a “first column” one. In this case, the proposed fusion formalism for the underlying lattice 
models [14,22] requires the following fusion rule for consistency: 

/C r> i x /Ci )S = JC r>s (r , s £ 2Z + — 1). (5.1) 

This is certainly consistent with the corresponding character product rule (3.15) and we have verified it explicitly, 
using the Nahm-Gaberdiel-Kausch algorithm, in many cases (see below). The evidence is, in our opinion, 
sufficient to conjecture that (5.1) holds in complete generality. Whilst this accords with the proposed lattice 
fusion calculations, we view the result as also confirming, indirectly, that we have made the correct abstract 
definition for Kac modules. 

We remark that when confirming the fusion product (5.1), the case in which /C rjS is a corner type module, 
hence is semisimple, is the most computationally intensive. To explicitly verify that each composition factor 
splits off as a direct sum, thereby forming the required collection of islands (as indicated in Figure 6), we 
must compute to the depth given by the maximal difference between the conformal weights of consecutive 
composition factors (when they are ordered by their conformal weight). This quickly becomes infeasible with 
our implementation as r and s grow, so the direct evidence for corner type modules is somewhat less compelling 
than for the other cases. 

We also note here that the computational complexity of the fusion algorithm means that we were only able 
to successfully confirm (5.1) when the required depth was at most 4. However, as p or pf increases, the labels r 
and s requiring a given depth calculation tend to increase leading to an overall steady decrease in the feasible 
depths due to the increasing complexity of the singular vectors for high r and s. For (p,p r ) = (2, 8) and (3, 5), 
we were therefore limited to depths at most 2. 

Including parity in the fusion rule (5.1) is easy: In each case, explicit computation confirms that 

IC+i x /C+ s = IC+ S (r, s € 2Z + — 1). (5.2) 

More generally, if the parities of JC r ,i and K. y s coincide (differ), then that of K. r<s will be even (odd). This 
observation does not appear to have a simple explanation in terms of the fusion algorithm, although it is in 
accord with the well known principle of conservation of fermion numbers. We will show in [12] that it follows 
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readily from a fermionic version of the Verlinde formula. From the lattice, it follows from considerations of the 
parity of the system size [14]. 

To illustrate some of the simpler issues that arise with fusion computations of the form (5.1), we consider 
the example £ 3,1 x K-\^, for (p,p r ) = (2,4) (c = 0). The Verlinde formula tells us that the character is that of 
/C 3,3 which means that the fusion product has two composition factors, £ 1/2 and £ 3 . However, there are three 
inequivalent structural possibilities: 


1. 
2 ' 


3: 




1 






(5.3) 


The first possibility is £ 3 , 3 , the expected result, the second is its contragredient dual, and the third is the direct 
sum £ 1/2 ® £ 3 - Constructing K, 3,1 x /Cy 3 to depth 0 leads to a one-dimensional truncated space. Comparing 
with the possible structures, we see that this is consistent with /C 3j 3 where the £3 factor is descended from C 1 / 2 , 
hence does not appear at depth 0. Moreover, the depth 0 truncations of the other two structural possibilities 
are two-dimensional. For this reason, a depth 0 calculation alone is sufficient to confirm that /C 3,1 x /C 1,3 = /C 3 , 3 . 

A slightly more complicated example is the (p,p r ) = (1,3) (c = —5/2) fusion product £ 3,1 x £ 1 , 5 . The 
character is that of £ 3 , 5 , so the composition factors of the product are £ 1 / 2 , £ 5/2 and £ 4 - Depth 0 computations 
reveal a two-dimensional truncated space with conformal weights 1/2 and 5/2. It follows from the general 
structure theory that £4 is descended from £ 5 / 2 , so it only remains to decide whether £ 1/2 splits off as a 
direct summand or whether it is generated from £ 5/2 through the action of the positive modes. This requires 
computing to depth 2 and the result indicates that £ 1/2 is not a direct summand. The structure is therefore 


(5.4) 


4 • 


confirming the fusion rule K, 3,1 x /C 1,5 = K, 3 , 5 . 

In the braided case, there may be further obstacles to overcome in completely identifying the structure of the 
fusion product. For example, the product /C 34 x /Cy 5 at (p,p') = (2,4) (c = 0) involves six composition factors: 
£ 0 , £ 1 / 2 , £ 3 / 2 , £ 3 , £5 and £ 21 / 2 - A depth 0 calculation indicates that all but £0 and £ 1/2 are descendants and 
a depth \ calculation shows that G 1/2 maps the weight \ vector to that of weight 0. The (partial) structure of 
the fusion product is thus 


(5.5) 



where it only remains to determine if there are upwards-pointing arrows emanating from the three lowest nodes. 
The N = 1 version of the Projection Lemma of [54, Lem. 5.1] rules out any such arrow from the node labelled by 
£- and from the nodes labelled by 3 and 5 to that labelled by 0. There are thus only two possible arrows: those 
from 3 or 5 to |. The presence of these arrows may be ascertained as in [31], see also [22, Sec. 4.2.2], If they 
are absent, then the fusion product would possess a singular vector of weight 3 or 5, respectively. By computing 
within the most general abstract module with structure (5.5) (an arbitrary extension of the higliest-weiglit 
module V 1 / 2 /V 15/2 by the highest-weight module /C 1.1 = V 0 /V 1 / 2 ), we can explicitly verify that such singular 
vectors do not exist. This is a relatively efficient calculation for a computer; in particular, we do not need to 
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invoke the Nahm-Gaberdiel-Kausch algorithm. In this way, we arrive at the structure 

o 

3 
2 

(5.6) 

5 


21 

2 


from which we conclude that £ 3,1 x /Ci ,5 = £ 3 , 5 . 



5.2. Fusing near the edge. We do not have enough data to make conjectures concerning the general fusion 
rules of the Kac modules (see the following section for some complicated examples). However, we have observed 
some patterns that seem to be followed when the Kac modules to be fused lie sufficiently close to the edges of 
the extended Kac table. Below, we will indicate what this means precisely. For brevity, we will refer to this 
situation as fusing “near the edge” (of the extended Kac table). 

Recall the following specialisations of the Kac character rules (3.14): 

r+r— 1 

ch[/C r> i] Elch[/C r / jS /] = ch[/C r » jS /], (5.7a) 

r"=\r— r'|+l 
s+s 7 —1 

ch[/Ci jS ] E3 ch[/C r / jS /] = ch[/C r / >s //]. (5.7b) 

s // = |s—s / 1 -|-1 

As usual, primed sums indicate that the index increases in steps of two. If the Kac modules /C r ", s ' (/C r ', s ") 
that appear in the above decompositions all satisfy either r" ^ p or s' Sj p' (r' St p or s" ^ p'), then we 
conjecture that the fusion rule corresponding to Equation (5.7a) (Equation (5.7b)) may be determined through 
the following procedure: 


1) Write down a list of all the Kac modules /C r ", s ' (AC r ', s ") from the decomposition (3.14) in order of increasing 
r" («"). 

2) Starting from the smallest value of r" ( s"), check whether there exists a /C p ", s ' (/C r < itT ") in the list which is the 
reflection of AC r //, s ' (/C r ',s") about the next boundary. This means that p" (a") must satisfy 0 < p" — r" < 2 p 
(0 < 0 " — s" < 2 p') and Alb ( r // +p ,/), s / (/C r ,, 1 ( S "+ CT ")) must be of boundary or corner type. 8 

3) If there does, then replace )C r » )S > and /C p ", s ' (/C r q s " and lC r p a ") in the list by the staggered module 
^’i( p // + r") s' (^r' HV'+*"))' -^ny logarithmic coupling must be determined through other means. 

4) Repeat with AC r ", s ' (/C r q s "), where r" (s") is the next-highest value. Once all values are exhausted, the list 
consists of the direct summands of the fusion product. 


Similar conjectures were made for certain classes of Virasoro Kac modules in [26,27,29,30]. We have checked 
that this procedure gives results that are consistent (up to the values of any logarithmic couplings) with our 
explicit fusion computations. Appendix C lists, for each central charge considered (except c = | for reasons 
that are explained in Section 5.5), those computations for which these checks have been performed. 

We illustrate this procedure with a few examples. First, take (p,p') = (1,3) (c = —|) and consider £ 1,3 x /C 1,5. 
The list, for this product, is /C 1 , 3 , /Cqs, K, 1 , 7 , according to (5.7b). Since /C 1,3 is of corner type, it has no reflection. 
However, 1C 1,5 reflects onto 1C 1.7 about the (Ramond) corner type module 1C 1,6, so they are replaced by <Sq’g. 
Computing the logarithmic coupling using Nahm-Gaberdiel-Kausch fusion or (B.5) then gives 


ICi, 3 x Kx, 5 = lCi, 3 © Si°6(-2) (c = -|). 


(5.8) 


Similar arguments for (p,p') = (2,4) result in 


1C 3,1 X /C3,3 = 5:2,’3(5) ©AC5,3, ACi ; 5 x /Ci, 5 


= 5° 4 (—4) ® ^(-l) © /Cp 9 (c = 0). 


(5.9) 


8 lt is worth mentioning here that this boundary or corner type module may well be from the Ramond sector. 
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We mention that in these examples, the original list of Kac modules did not by itself uniquely determine which 
list members are combined to form a staggered module. This is where it is important to start the above 
procedure with the smallest value of the appropriate Kac label. For example, with ACy 5 x AC 1 , 5 , we combined 
/C 14 with ACy 7 to correctly identify <S 4 4 (—4) as a direct summand, instead of combining ACy 7 with /Cyg. 

The parities of the modules obtained by fusing near the edge of the extended Kac table are easily determined. 
Assuming that the Kac modules being fused are both assigned an even parity, we find that the parities of the 
Kac modules in the ordered list constructed in step 1) above always alternate, starting (and therefore ending) 
with even parity. This is also consistent with lattice expectations [14]. As an example, take {p,p') = (2,4) 
(c = 0) and consider the fusion product ICf 3 x /C]~ 4 . This yields the list /C^, ^“ 2 ~ 4 > ^26 f rom which we deduce 
that the first and last Kac modules combine to form a staggered module. As ft. 2,2 = /i 2,6 and both /C 2 ,2 and 
AC 2,6 are highest-weight modules (see Figure 5), there is no logarithmic coupling to find and the final result is 

/C+ 3 x ^2)4 = <S°,4 + © ^ 2,4 (c = 0). (5.10) 

As in the previous section, this parity rule will be derived in [ 12 ] from a fermionic Verlinde formula. 

The (conjectural) procedure described above for fusing near the edge of the extended Kac table implies that 
the resulting fusion products may only admit Jordan blocks of rank at most 2 for the action of Lq. More precisely, 
it implies that these products always decompose as direct sums of Kac modules and staggered modules. We 
have observed such staggered modules in every model considered except for that with (p,p r ) = ( 1 , 1 ) (c = |). 
This exceptional case is discussed separately in Section 5.5 below. 

Finally, the procedure proposed above, for determining Kac fusion rules near the edge of the extended Kac 
table, does not require restricting to fusing with either AC r ,i or ACy s , although the Kac module labels (■ r",s ") 
appearing in the character product (3.14) may need to satisfy either r" ^ p or s" ^ p'. An example illustrating 
this is {p,p') = ( 2 , 4 ) (c = 0 ) with /C 2 ,2 x AC2,2- The corresponding list in this case is /Cy 1 , £3,1, /Cy3, £3,3 and 
an explicit calculation shows that the first two and last two members combine to form staggered modules as 
one would expect from the boundary reflection principle: 

IC 2 , 2 X /c 2 ,2 = S^d) ® $!#(§) (C = 0). (5.11) 

Further examples like this are common for larger p and p' and we have checked in several cases that the fusion 
decompositions do lead to staggered modules whenever two Kac modules in the list (not ordered in these 
examples) are related by reflection about a boundary, see Appendix C. However, the number of examples that 
we are able to fully analyse is not particularly large, explaining why we have not included these observations in 
the conjectured procedure. 


5.3. Fusing away from the edge. We first consider an interesting product, AC 1,3 x /C 2)2 at (p,p r ) = (1,3) 
(c = — f), which takes us slightly out of our comfort zone near the edge of the corresponding extended Kac 
table. The Kac module list given by (5.7b) is /C 2 , 2 , AC2,4 (since the character of AC2,o is formally 0, by (3.8)) 
and we note that ( r',s") = (2,4) is not near the edge — it fails to satisfy both r' < p = 1 and s" ^ p' = 3. 
However, K. 2,3 is of corner type, so one might expect that the result is the staggered module <S 23 (with some 
logarithmic coupling). It is, at first, surprising that the fusion product is actually found to be the staggered 
module <S 4 g(—2) of (5.8). However, there is no contradiction as the structures are identical: 



(5.12) 


The fact that G_ 1/2 annihilates the weight 0 vector of AC 2,4 means that its action on the preimage of this vector 
in <S 2 % will be proportional to the weight J vector of the AC 2.2 submodule. The fusion algorithm merely shows 
that the proportionality constant is not zero. In this case, the conjectured procedure of the previous section 
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does predict the correct answer. However, it indicates strongly that a theory of staggered extensions of Kac 
modules, rather than highest-weight modules, will be needed to properly understand the Kac fusion rules away 
from the edge. 

The next example is a much more structurally intricate fusion product for which the conjecture of the 
previous section fails spectacularly. Here, we consider K.2,2 x £2,4 at {p,p') = ( 2 , 4 ) (c = 0 ), the result of which 
is sufficiently complicated that we have not unravelled its full structure. First, the corresponding character rule 
suggests that the product should involve the following Kac modules (which take us well away from the edge of 
the extended Kac table): 



( 5 . 13 ) 


In particular, this implies that the composition factors (including multiplicities) of the fusion product are 2 £o, 
4 C1/2, £3/2, 2 £3, 2 £5 and £21/2- If the result were the direct sum of these Kac modules, then we would see 
five linearly independent eigenvectors when computing with the Nahm-Gaberdiel-Kausch algorithm to depth 0 , 
three of eigenvalue I and two of eigenvalue 0 . However, the special subspace of /C2,2 is only four-dimensional, 
so such a direct sum is ruled out. Indeed, explicit computation shows that the depth 0 truncation of the fusion 
product has conformal weights 0 and I, both with multiplicity 2 , and that Lq possesses two rank 2 Jordan 
blocks. 

Extending to depth I, we encounter a new feature: our first rank 3 Jordan block. The placement of the £3/2 
factor is now determined by looking at the staggered submodule generated by the copies of £0 and applying 
Proposition B .3 to rule out one of the two possibilities: 



( 5 . 14 ) 


We have confirmed this placement by computing to depth |. However, determining the placement of the £3 
factors using the fusion algorithm is too computationally demanding at present. We are therefore left with the 
question of where to place the composition factors 2 £3, 2 £5 and £21/2- The first thing to note is that £21/2 
can only appear as a descendant of both an £3 and an £.5. We are able to further restrict the possibilities by 
appealing to staggered module theory. 9 The Jordan block structure thereby narrows the structural possibilities 
down to just two (though this is before upwards-pointing arrows are considered): 



( 5 . 15 ) 


9 Hcre, we look at the constraints on the existence of the staggered subquotients of this module. Some of these are simple conse¬ 
quences of Proposition B.3, but more powerful non-existence theorems follow from the Neveu-Schwarz analogues of [54, Sec. 7]. 
Similar arguments (in the Virasoro case) may be found in [22, Sec. 4.2.2]. 
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However, we note that even if we were able to decide between these structures (by computing to greater depths, 
for example), then it is still not clear if there are additional parameters, generalising the logarithmic couplings 
of staggered modules, required to completely identify the isomorphism class of this indecomposable module. 

The existence of indecomposable Neveu-Sclrwarz modules like these is not unexpected given that similar 
modules have been investigated over the Virasoro algebra [ 29 ]. However, a full understanding of a single 
example of these modules is still missing (see [ 48 , Sec. 6 . 4 ] for a recent discussion), hence attempts to further 
explore these structures here would be misguided. 

Nevertheless, an obvious feature of /Cj.3 x /C2,2 at c = 0 is that the Kac modules ( 5 . 13 ) suggested by the 
character product are related by boundary reflections of different orientations. For example, /C+3 reflects onto 
/Ci,5 about the boundary type (Ramond) module /Cq4, but it also reflects onto £3,3 about the boundary module 
AC2,3. With this observation, it is not surprising that the procedure conjectured in Section 5.2 does not suffice. 
It would be interesting to further investigate whether this procedure gives correct fusion results whenever this 
multiple orientations phenomenon is absent. An issue to address here is what happens when the (expected) 
structure of the fusion product includes non-cyclic Kac modules, meaning that the module is not generated by 
a single subsingular vector. Non-cyclicity is generic for Kac modules (see Figure 6), but current computational 
limits mean that our fusion calculations do not shed any light on this issue. 

5 . 4 . Fusion subrings. As a (finitely generated) fusion ring is commutative and associative, it may be presented 
as a quotient of the polynomial ring over Z in which the indeterminates are the generators. In this section, we 
derive conjectures for such presentations of the fusion subrings generated by either /C3,1 or /Cq3, generalising 
similar conjectures for Virasoro fusion subrings made in [66]. 

We begin by discussing polynomial presentations for the ring of Kac characters equipped with the Verlinde 
product of Section 3 . 2 . In this case, the character product ( 3 . 14 ) exhibits an obvious sl( 2 ) symmetry that 
suggests a presentation in terms of Chebyshev polynomials. Indeed, it is easy to check that the subrings 
generated by ch[/C3,i] and chf/Ci^] are isomorphic to Z[x\ and Z[y\, respectively, and that explicit isomorphisms 
are given by 

ch[/C2i+i,i] <-> U2i(\^/x + l), ch [/Cl ,2j+l] ^ U'2j (5 \J y + l) > ( 5 . 16 ) 

where U n denotes the n-th Chebyshev polynomial of the second kind. Despite the square roots, the images of 
the Kac characters are polynomials of degrees i and j, in x and y , respectively. 

We remark that the subring generated by ch[/C3,i] and chf/Ci^] is obviously isomorphic to Z [x,y], but that 
this does not contain the Neveu-Sclrwarz Kac characters with even indices. To rectify this, one may introduce 
a new generator z corresponding to ch[/C2.2] and show that the full Kac character ring admits the following 
presentation: 

(ch[£3.1] 1 ch [K.2,2] 1 ch[£1,3]) a ( 5 - 17 ) 

Before we can lift ( 5 . 16 ) to presentations of the fusion subrings generated by either /C3,1 or /Cq3, there are 
a few matters to address. First, we mention that p = 1 implies that the fusion subring generated by /C3,1 is 
semisimple (see Section 5.5 below). In particular, it is isomorphic to the character ring generated by ch[/C3,i]. 
The story is identical for p' = 1 and /Cq3. Second, the case p = 2 also requires a separate treatment because 
fusing two (non-trivial) elements of (£3,1), such as /C3,1 x £3,1, never gives /C3,1 as a direct summand. Again, 
the story is identical for pf = 2 and /C 1,3. For simplicity, we will therefore assume that p ^ 3 or pf Js 3 when 
analysing the fusion subring generated by /C3,1 or /C 1,3, respectively. 10 

From here on, we will restrict to the fusion subring generated by /C3+ as the results for the subring generated 
by /Ci,3 then follow by interchanging labels and swapping p for p'. Consider the n-fold iterated fusion product 
/C 3 " = /C3+ x ••• x /C3+. The (conjectured) procedure of Section 5.2 states that for n ^ this iterated 

fusion product will decompose as a direct sum of Kac modules /C21+1.1, with i £ { 0 , 1 ,..., }• However, 

when n = > the fusion product will have a staggered module as a direct summand. Unfortunately, we have 

not computed any fusion rules in which one of the modules being fused is staggered. 

10 The mild pathology exhibited when p or '[/ is equal to 2 is a manifestation of the fact that we are ignoring the Ramond sector 
fusion. Nevertheless, the results of this section can also be applied in these cases with minor modifications. 
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If p is sufficiently large, however, the fusion rules of the staggered modules may be derived using associativity. 
For example, if p ^ 5 is odd, then 

^03,1 x K-p ,l = fcp, 1 ©<Sp’i 

implies that 

/C 34 x = (/C 34 x ^-3,1 x /Cp,l) 0 (^-3,1 x fcp, l) = /Cp, 1 0 (/C. 5,1 x /Cp,l) = 2 /Cp,i © <Sp’i ©<Sp’i, 
where we ignore any logarithmic couplings. 11 A similar calculation shows that if p ^ 6 is even, then 


^3,! x Sp’° = 2 ^ 0 Sp’i- 


1,0 


-3,0 


(5.18) 

(5.19) 

(5.20) 


One can also obtain the fusion rules involving other Kac modules and those of staggered modules with one 
another (for sufficiently large p). Calculations of this nature lead us to conjecture that the indecomposable 
modules appearing in the fusion subring generated by £34 are precisely as follows: 

(£3,1) = span z |/C 2 i+i,i,/C m p > i,5^ 0 1 : i = 0 , 1 ,..., ; fc,rnGZ + ; 0 < a < p; kp - a = mp = 1 mod 2 j. 

(5.21) 

Here, we again ignore any logarithmic couplings. Note that if p is even, the /C m p,i do not belong to this subring. 

Combining the (conjectured) results of Section 5.2 with associativity leads to a general conjecture for the 
fusion rules of the (unital) subring generated by /C 34 . This is conveniently encoded in the following polynomial 
ring structure, generalising that presented above for the corresponding character subring. In particular, we lift 
(5.16) from characters to the Kac modules in (5.21): 

^V 2 z+i,i ^ ^ U 2 i {^~2 V 4: © i() (/ — 0 , 1 ,..., [_■ 2 J), ^dmp,i aa U m p—i (^\/x T1) (tti G Z.^., vnp — 1 mod 2). 

(5.22) 

To obtain the map for the staggered modules S^ p \ in (5.21), we recall that (5.16) maps their characters to 
ch[«S^°i] = ch[/C k p-a,i] + ch[/C fcp+a ,i] 

AA Ukp-l-a(^Vx + l) + Ukp-l+a(^Vx + l) = 2 T a (^y/x + l)Ukp-l x + l), (5.23) 


where T n denotes the n-th Chebyshev polynomial of the first kind. Lifting this to modules, we arrive at 

Sk p °i aa 2T a (±->/x+ l)Uk P -i(\y/x + l) {k G Z + , 0 < a < p, kp - a = 1 mod 2). (5.24) 

As the degrees of the images in (5.22) and (5.24) are i, \{mp — 1) and \{kp — 1 + a), respectively, and as the 
coefficient of the highest degree term of each image is 1, we see that these maps define a Z-module isomorphism 

from the fusion subring (£ 3 , 1 ) to Z[x]. Our conjecture for the fusion rules is that this is actually a ring 

isomorphism. Similarly, we conjecture that the map from (/Ci, 3 ) to Z[y], defined by 

fci,2j+i aa U2j(\^/y +1) 0’ = 0, l,..., 

A-i.np 7 AA U np '-i(\sjy + l) {n G Z + , np' = 1 mod 2), (5.25) 

s Up' 0 ( 2 ~ $>=o) T b(\\/y+ l)Ue p '-i(^y/y + l) (e G Z + , 0 < b < p', ip' - b = 1 mod 2), 

is a ring isomorphism. 

To illustrate these conjectures, we consider how (5.25) allows us to compute the fusion rule /C 1,3 x /C 1,3 when 
C PiP') = (2> 4): 

/Ci , 3 x /Ci, 3 aa yy = 1 + (y 2 - 1 ) aa /C M 0 (c = 0 ). (5.26) 

The result agrees with (4.19), up to the logarithmic coupling of the staggered module. It is very easy now 
to predict more complicated fusion rules, as long as the modules being fused belong to either of the subrings 
considered above. For example, (5.25) predicts the (p,p') = (2,4) fusion rule 

<S°4 x $$ ( y 2 ^ l)(z/ 3 ^ 2y 2 - y + 2) 

= 2 (y 2 - 1) + 2 (y 4 - 2 y 3 - 2 y 1 + 2 y + 1) + ( y 5 - 4 y A + 2 y 3 + 6 y 2 - 3 y - 2) 

llr rhere are further questions to address here concerning whether the coupling /3^’J obtained, for example, from /C 3 ,i x /Cp,i is the 
same as that obtained from £ 5,1 x /C p> i. The evidence that we have collated to date indicates that the answers are always yes. 
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o 2 S 0 ^} © 2 5° g 1 © (c = 0). (5.27) 

We remark that these conjectures may be incorporated into the procedure described in Section 5.2 for fusion 
near the edge of the extended Kac table. If either of the modules to be fused is staggered, say, then one 
replaces it by /C r - a , s -b ® fcr+a,s+b when applying the character decomposition (3.14) in the first step. 


5.5. An exceptional case: c — In the exceptional case where {p,p') = (1,1), every Kac module has the 
islands structure, hence they are all semisimple. The fusion rules therefore reduce to those for the simple Kac 
modules and these are exhausted, up to isomorphism, by the /Ci jS (in particular, JC r> i = /Ci jr ). According to 
the conjectured prescription in Section 5.2, the fusion products of these simple modules always decompose into 
a direct sum of Kac modules K,\. s " because corner type modules do not have reflections about boundaries. It 
follows that, in this case, the general character rules (3.14) lift to the genuine fusion rules 

r+r 7 —1 s+s 7 —1 

/+« x K r ,, s , = /C r " jS " (c = §), (5.28) 

r"=\r— r'|+l s" = |s—s'| + l 


where the primes on the direct sum symbols indicate that the summation variables increase in steps of two. 
The fusion coefficients corresponding to this rule were previously derived in [67]. 

We remark that one can derive this result from the isomorphism /C rjS = (specific to this model), the 
conjectured fusion rule (5.1), and the definition of corner type Kac modules. Indeed, we have 

r+s—1 r+s 7 —1 r+s—1 r+s 7 —1 

/C r , s x JC r f }S / = ^ (J) A+x ^ /Ci jCr ^ = (J) ICp.a 

p=|r-s| + l <7=\r' —s' | + 1 p=|r-s| + l cr=|r'-s'| + l 

r+s—1 r'-\-s' — 1 p+cr—1 r+r 7 —1 s+s 7 — 1 r"-\-s" — 1 

= ® © ® /Cl.<r' = © © © /Cl, v" 

p=|r—s| + l a=\r'— s'| + l cr' = |p— cr| + l r" = |r—r'| + l s'^ls—s'| + l cr" = \r" — s" | + 1 

r+r 7 —1 s+s 7 —1 

= /C r ", s ". (5.29) 

r"— |?—r'|+l s // = |s—s'|+l 


Of course, semisimplicity also proves that the character product (3.15) lifts to the fusion rule (5.1), when c= |. 
It follows, in particular, that the Kac modules form a closed fusion ring, without the need to introduce any 
staggered modules. This does not mean, however, that staggered modules do not exist at c = |. For example, 
one can construct staggered self-extensions of every (simple) Neveu-Schwarz Kac module /Ci jS , with s / 1, by 
following the arguments of [54, Prop 7.5]. 


6. Discussion and Outlook 

The results presented in this article demonstrate that the paradigm of fusion product computations using the 
Nalnn-Gaberdiel-Kausch algorithm is as successful for the Neveu-Schwarz algebra as for the Virasoro algebra. 
In particular, our computations have allowed us to formulate several conjectures for Kac module fusion rules in 
the logarithmic N = 1 superconformal minimal models. We have seen that staggered Neveu-Schwarz modules, 
upon which Lq has rank 2 Jordan blocks, are readily encountered and that more complicated indecomposable 
modules can be generated on which Lq has rank 3 blocks. These results provide, among other things, strong 
evidence for the conjectures made in [14] from numerical lattice-theoretic studies. On the other hand, the 
reader may have noticed that our results and conjectures bear a striking resemblance to their counterparts for 
the Virasoro algebra, see [22,27,29,30]. This accords with expectations, as there are many instances in which 
the representation theories of the Neveu-Schwarz and Virasoro algebras mirror one another. 

From this point of view, the key advance made in this paper is not so much the explicit results themselves, 
but rather the development of the formalism that led to them. This should be understood in the context of the 
extension of this formalism, to include the Rarnond algebra, that will appear in [12]. In this case, the standard 
module development that led to the logarithmic Verlinde formula will need to be generalised to accommodate 
supercharacters, ultimately resulting in a fermionic Verlinde formula generalising that of [68]. More tellingly, 
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the Nahm-Gaberdiel-Kausch fusion algorithm will require significant modifications because the Rarnond sector, 
consisting of twisted modules, forces one to deal with non-integer field-theoretic monodromies. The required 
changes have been partially addressed in [36], but this work requires refinement (and simplification) before the 
twisted version of this algorithm can be efficiently coded and utilised. 

Let us recall that the representation theory of the Rarnond algebra can be significantly more involved than 
that of the Virasoro or Neveu-Schwarz algebras. In particular, the relation Gq = Lq — c/24 means that Go may 
act nilpotently, hence possibly non-semisimply, on weight spaces. Such a non-semisimple Go-action is observed, 
for example, on the Rarnond Verma module of highest-weight c/24. Because it seems that the majority of the 
difficulties in the Rarnond sector may be traced back to this observation, and because nilpotent actions are 
the norm for the odd generators of the N > 1 superconformal algebras, it is reasonable to regard the N = 1 
Rarnond sector as a non-trivial, but accessible, toy model for the representation theory of these important 
superalgebras. This gives a separate motivation for studying logarithmic N = 1 models: They give us an idea of 
what one can expect in the N > 1 case. As logarithmic behaviour is now recognised to be generic (rather than 
pathological) for conformal field theories, it makes sense to look for logarithmic structures when investigating 
poorly understood representation theories. 12 

We mention that such superconformal algebra studies fit rather naturally into a bigger research programme 
that is currently being actively pursued (see [70] for further details). Here, the point is that the superconformal 
algebras may be naturally constructed as quantum hamiltonian reductions of certain affine Lie superalgebras [71]. 
For the N = 1 superconformal algebras, the corresponding affine Lie superalgebra is osp(l|2). The conformal 
field theories corresponding to these affine superalgebras are expected to be logarithmic, in general, and the 
relatively tractable affine symmetry suggests that these theories will be very valuable for understanding general 
logarithmic behaviour. Unfortunately, little is known about these models (but see [42,72-76]). A more thorough 
understanding of these affine theories is therefore warranted and we expect that the structural features of the 
logarithmic superconformal models will aid in this undertaking (and vice-versa). 

This work also relates to affine models through the celebrated coset construction of Goddard, Kent and 
Olive [77,78]. It is well known that the simple modules of the N = 1 minimal models arise when considering 
admissible highest-weight modules [79,80] of the sl(2) components of the coset. However, as proposed in [81], 
one expects to similarly realise N = 1 Kac modules, in particular, by extending these considerations to the 
analogues of Kac modules over sl(2). As this approach is also applicable to the infinite hierarchy of extended 
minimal models [82-87], this will give important insight into their logarithmic counterparts. These ideas are 
also currently being actively pursued. 

Finally, it would be of interest to explore whether our results admit a W-extended picture, generalising 
the situation [25,33,88-95] for the Virasoro logarithmic minimal models. From the lattice [96-98], this would 
amount to identifying a W-extended vacuum boundary condition for the Neveu-Schwarz algebra and devising 
the appropriate lattice implementation of fusion of the W-extended representations. Presumably, this would 
give rise to W-extended Neveu-Schwarz representations whose characters could then be compared with the 
recent results of Adamovic and Milas [19,99,100]. 
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Appendix A. Fusion and the Nahm-Gaberdiel-Kausch Algorithm 

In this appendix, we review an algorithmic approach to fusion that was proposed by Nahm [50] and then 
significantly generalised by Gaberdiel and Kausch [28]. This approach aims to construct the fusion product 

12 We mention an example of where this looking has paid off: The longstanding issue of negative “fusion coefficients” for fractional 
level sl(2) models, see [43,69]. 
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of two (vertex operator algebra) modules by realising it as a vector space, in fact as a quotient of the vector 
space tensor product of the modules being fused, upon which the chiral algebra acts through explicitly given 
coproduct formulae. These formulae are deduced from fairly straight-forward manipulations involving operator 
product expansions [52,101] and the aforementioned quotient corresponds to imposing the condition that two 
seemingly different means of arriving at these coproduct formulae actually give identical results. As we will see, 
this imposition follows from the mutual locality requirement for conformal fields. 

The algorithm that has come to be known as Nahm-Gaberdiel-Kausclr fusion then observes that this tensor 
product realisation can be consistently truncated so as to explicitly construct only a quotient, preferably a finite¬ 
dimensional one, of the fusion product vector space. It is often the case that one can completely identify a given 
fusion product by analysing the coproduct action on a sufficiently large truncation. This allows one to perform 
fusion calculations explicitly with a computer algebra system. There are several advantages to this approach 
to fusion over traditional approaches, the most important of which is that it facilitates the exploration of the 
new classes of representations that fusion produces. For example, the result of fusing two simple modules need 
not be a direct sum of simple modules in general; in fact, the fusion product need not even be highest-weight. 
The main disadvantage is that the algorithm amounts to a computationally intensive brute force construction 
and is therefore not well suited for general fusion calculations or theoretical studies. There is also the issue 
of determining when one has determined all of the so-called spurious states. The algorithm itself does not 
guarantee termination here, so it is important to be able to check this independently. In this article, we employ 
a Verlinde formula for this purpose; in some cases, one can instead use correlation function computations to 
confirm the result a posteriori. 

Readers who are already familiar with the Nahm-Gaberdiel-Kausclr fusion algorithm, as introduced in [28], 
will find little that is new here, though they may wish to skim what follows in order to familiarise themselves with 
our notation and nomenclature. The Neveu-Schwarz supersymmetry does lead to the appearance of additional 
signs, as one expects with a graded tensor product, but the ideas and implementation follow the bosonic case. 
On the other hand, fusing with a Ranrond module requires significant modifications to the algorithm. We shall 
defer a detailed discussion of these modifications to [12]. 


A.l. Coproduct Formulae. Here, we review the derivation, following [52], of the coproduct formulae that 
define the action of the chiral algebra on the fusion product. We provide almost all of the details because the 
generalisation to the twisted sector, necessary for Ranrond fusion, is significantly more involved. In fact, our 
treatment of the twisted case, which will be reported in [12], simplifies the known formulae [36] considerably. 
We mention that the derivation for negative modes is only valid in the limit when one of the insertion points is 
sent to zero. This is not important for practical purposes, but is expected to be relevant to demonstrating that 
these fusion coproducts define a tensor structure on appropriate categories of vertex operator algebra modules. 
Let ( 2 ) be chiral fields of conformal weight whose mode expansions take the form 

£(%)= ^ S^z~ n ~ hW . (A.l) 

ngZ—ftW) 

For the purposes of this article, the relevant fields are T(z), with h = 2, and G(z), with h = |. In what follows, 
we will omit the index j labelling the chiral field for simplicity. We also introduce arbitrary fields ipi(wi ) and 
^ 2 (^ 2 ) that are local with respect to 5 ( 2 ): 

S(z)ipi(wi) = piipi(wi)S(z). (A.2) 

Here, € C\{0} is the mutual locality index of S with tpi (this index also depends on S; the notation likewise 
keeps this implicit for simplicity). 

We deduce coproduct formulae for fusion by determining the natural action of the modes S n on the (radially 
ordered) products ipi (u>i) { w 'i) ■ This will define an action of the S n on the tensor product of the corresponding 

states ipi <g> if) 2 , the latter being interpreted (after quotienting) as a state in the fusion product. The starting 
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point of the computation is the contour integral 




where the contour F encloses 0, uq and w- 2 , fl is the vacuum, and <f> is an arbitrary spectator state. We remark 
that cf) may even depend on other insertion points, noting that because of radial ordering, such insertion points 
are not enclosed by T. 

We will assume that the fields ipi{ w l) and ^ 2 (^ 2 ) correspond to untwisted representations of the chiral 
algebra, referring to [12,36] for the twisted case. In other words, their operator product expansions with S(z) 
are characterised by modes S m with m € Z — h: 

S(z)ipi(wi)= ^2 - Wi)~ m ~ h . (A.4) 

Inserting these operator product expansions into (A.3), we see that there are no branch cuts in the integrand. 
We may therefore split the contour into three simple contours around each of the (potential) singularities w±, 
W2 and 0. The residue at w\ is computed by substituting the operator product expansion of S(z) and 

<f {(/>,S(z)4> 1 (wi)'lp2(w2)tyz n+h ~ 1 ^ :: T = y 2 (</>, (Smi>l){wi)i>2{w2)fl) <f (z - W 1 )~ m ~ h Z ^* 1-1 ^7 

Jwi 2th , J w , 2tti 


= E 


m=-h -\-1 


n + h — l 
m + h — 1 


Wi m (<j>,(S m 'ipi)(wi)ip2('W2)ty- (A.5) 


The residue at W 2 is similarly computed, though we must first apply (A.2) before we are able to substitute the 
operator product expansion S[z)ip 2 (w 2 )- 


((j),S(z)'t/j 1 (w 1 )lp2{’W2)^)z n+h 

J Wo 


_-j dz 


n + h—l 
in + h — 1 


u>2 m (</>, 4>i{wi){S m il)2)(w2)V). (A.6) 


Note that (A.6) follows readily from (A.5) upon swapping wi with n >2 and then making the replacements 

S m ipi -+ V’l and i >2 -+ liiSmfo. 

If n ^ —h + 1, then the integrand of (A.3) has no pole at z = 0 and the coproduct formula for S n is just 

a wl , W3 (s„)= (»>-/»+ 1 ). ( a . 7 ) 

Here, we have extracted the action of S m on the fields in the correlator so that, for example, 

(</>, {S m i/Ji)(wi)il>2(w2)ty —* (S m ® 1). (A.8) 

It should be clear now why we have insisted on the spectator state <f> even though it plays no role whatsoever in 
the analysis. If n ^ —h, then there is a pole at z = 0 and we can evaluate the corresponding residue by using 
either operator product expansion: 

<j>{(l),S{z)%l) 1 (w l )'ip 2 {w 2 )^)z n+h ^ 1 = S'] l(z-wi)“ ro_ V +i, “ 1 ^ : (S ra ®l) (A.9a) 

= 1+ V £(z-W 2 )- m - h z n+h - 1 p T (l®S m ). (A.9b) 

g-rjf h JO 

We remark that 2 is supposed to make a small circle around 0, while inserting the appropriate operator product 
expansion requires that 2 be close to either w\ or uj- 2 . We conclude that (A.9a) is only valid when w\ is close to 
0 whereas (A.9b) is only valid when u >2 is close to 0. We will later send w\ or W 2 to 0 to arrive at a simplified 
form for the fusion coproducts. 

Evaluating the integral in (A.9a), we see that the 2 = 0 contribution to the fusion coproduct A Wl>W2 (S n ), 
with n ^ —h, takes the form 


—m — h 
—n — h 


(-wO^iSm ® 1 ) 
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™ A ' m =—oo a ' 


m-—h-\-1 

Here, we have used the binomial coefficient identity 


—in — ft 
—n — ft 


= (- 1 ) 


-n-h I'm-n - 1 
—n — ft 


= (-l) 


-n- h fm-n- 1 
to + ft — 1 


= (-l) 


m— n—1 f n + h 1 

TO + ft — 1 


(A.ll) 


which is valid when to ^ —ft + 1 and n SC —ft. It is clear that the first sum on the right-hand side of (A.10) 
precisely cancels the contribution (A.5) from z = w i, when n ^ —ft. Similarly, starting from the evaluation of 
(A.9b) shows that the contribution (A.6) from z = w 2 is cancelled. 

In this way, we arrive at two formulae for the fusion coproduct of S„ with n ^ — ft: 

AV W 2 (S n )= (~^Zh) { - Wl)n ~ m{Sm ® 1) + ^ £ (rnt£lV’ rm(1 ® 5m) (A ' 12a) 

m=— oo A ™ — -h- i—i a ' 


h 


a $U(s n )= e (^ +h h %r m (s m ®i)+i*i E (A - 12b) 

m=-h-\-l ' ' m= — oo a ' 

Of course, these formulae should coincide in some sense. However, they are only supposed to be valid when w\ 
or W 2 is small, respectively, so we make the substitutions w 1 = 0 and W 2 = —w in the first and w\ = w and 
W 2 = 0 in the second to obtain 


OO / 

AW_ w (S n ) = (S n ®l) + vi E 

rr> — —h -LI A 

OO 

A $,($„)= E 


m=-h-\-1 

n + h—l 

TO + ft. — 1 


n + h—l 
to + ft — 1 


(-«y- ro (i ® s m ) 


w n ~ m (S m ®l) + Lii(l®S n ). 


(A.13a) 


(A.13b) 


m=-h-\-1 

We remark that we can set w\ = 0 or W 2 = 0 in the first or second coproduct formula above because the 
binomial coefficient (Z_™Zh) vanishes if to > n. We would not have been able to make these simplifications if 
the residue at 2 = 0 did not partially cancel the contribution from z = w\ or 2 = W 2 - 

To summarise, the fusion coproduct formulae (for untwisted representations) may be expressed as 

AoE(E) = {Sn 0 1) + hi E ( H + _ 1 ) ® Sm) (n>-ft + l), (A.14a) 


m=-h-\-1 
00 


A ( o,L(S-n) = ( S-n ® 1) + Ml E ( W 1 1 - (-l)” n+,l “ 1 ^- n - m (l ® 5 ro ) (n > ft) 

m=-h+l A / 


(A.14b) 


or as 


A$(S») = E 

m=—h-\-l 
oo 

Al 2 ) 0 (S_„)= E 

m=—h-\-l 


n + h—l 
to + ft. — 1 

to + n — 1 
to. + ft — 1 


w n - m (5 m 0 1) + M i(l®S ri ) 


(n ^ — ft + 1), (A.14c) 


(- 1 ) 


m+/i—1 —n—m 


w 11 m (S' m 0 1) + Hi{l 0 S-n) (n ^ ft). 


These formulae are related by translation: 

A£l w (S n ) = A^ 0 (e™ L -iS n e-™ L -i). 

Using [L— 1 , 5 n ] = (—ft + 1 — n)S n -\, we derive inductively the following formulae: 


LtiSn = E ( ■ ) + !-«)•■■ (-ft + j - n)S n -jL l _E 


t?L_ 


1=0 

Consequently, one has 

OO 

AgE(-S-n) = E 

1=0 


= E 

1=0 


—ft + j — n 

j 


(A.14d) 


(A.15) 


S n -je wL ~ 1 . 


(A.16) 


n — h + j 

j 


'<>„(«-.->) = E (“_ A”-”aL 2 ,US-J (A. 17 ) 

m—n ' ' 
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which leads to the three master equations for untwisted fusion: 


A(S„) = E {n>-h + 1) 

m=-h +1 ' ' 

(A.18a) 

A (S- n )= £ (n>h), 

ra= —/i+l ' 

(A.18b) 

S_„ <8 1 = ^ A(5_ m ) + , l( -l)-^ £ ( m ? t-7 1 W n-m ( 1 ® 5 ro ) (n > h). 

m=n ^ ' m=— h+1 ' ' 

(A.18c) 


/ O') 

Here, we let A denote A^, 0 for brevity. In practice, such as when performing the explicit computations reported 
in this work, we would set w to 1 to further simplify these formulae. However, this masks the natural grading of 
these formulae by conformal weight in the same way that choosing insertion points masks the conformal grading 
of operator product expansions. 

We remark that the interpretation of fusion as a quotient of the (vector space) tensor product is captured in 
Equation (A. 15). This formula is actually a requirement imposed by the locality of operator product expansions 
on the coproduct formulae. Thus, we may define the fusion product (as a vector space) of two chiral algebra 
modules M and N to be the quotient 


M x N = 


M ®c N 


((A£E($0 - A^ > o(e tui - 1 S„ e - wi - 1 )) (M ®c IV)) ’ 


(A.19) 


where the ideal is the sum of the images for all chiral modes S n (and all w). Of course, this is the same as 
the intersection of the corresponding kernels. The point is that M x N is defined to be the largest quotient of 
Ad ®c N upon which the coproduct actions coincide. We expect that this can be interpreted as a universality 
property for fusion as is imposed in [102]. 

Finally, we remark that the cancellation between the contributions to the integral (A.3) from the residues 
at z = 0 and z = w\ or z == w? is best understood analytically in terms of the regularity of the integrand at 
infinity. Explicitly, if we use the operator product expansion S(^z)^i(wi), then the sum of the contributions is, 
for fixed m ^ — h + 1 and n ^ —h, proportional to the integral 




-m-h^n+h—l 

2ir\ 


\ —m—n n -\-h—1 X ( \ — m—n 

W i) z n+n L — = -(b(z-w i) z 


—m-h^n+h -1 

2m 


= 0(1 - Wl y)- m - h y™-”-l 


= 0, 


(A.20) 


because m — n — 1 ^ — h+l + h— 1 = 0. In the above derivation, we have chosen to instead derive the cancellation 
of these contributions algebraically, using binomial coefficient identities, because this method readily generalises 
to the twisted sector [12], whereas the contour analysis becomes rather more subtle there due to the presence 
of branch cuts. 


A.2. The Nahm-Gaberdiel-Kausch Fusion Algorithm. One important insight [28,50] into the definition 
(A.19) of the fusion product is that it admits consistent truncations which are easier to construct explicitly. 
For this, we consider subalgebras it of the universal enveloping algebra of the chiral algebra. One important 
example is that generated by all the chiral modes with index not greater than minus their conformal weight: 

H ss = (s® : n < -h&y (A.21) 

For the Neveu-Schwarz algebra, this subalgebra is generated by the L n , with n < —2, and the Gj, with j < — 
The second example, actually a family of examples labelled by d £ Z, that we shall employ is that spanned by 
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monomials of weight greater than d in the chiral modes: 

il d = span {S^SW ■ ■ ■ Sfc) : r G Z^ 0 , «i + n 2 + • • • + n r < -d}. (A.22) 

The integer d will be referred to as the depth. 

The ability to consistently truncate the fusion product amounts to the following claim [28,50] for the chiral 
algebra modules M and N: 

M ^ AT 1VJ JV 

— (fss Tl\,f\ (Id! A T\ ’ (A.23) 


M x N 


il d (M x N) ~ ii ss {M) ^ ii d {N)' 

The hrst factor on the right-hand side defines the special subspace M ss of M and the second factor defines 
the depth d subspace N d of N (even though both are defined as quotients). The claim is therefore succinctly 
expressed as the inclusion (M x N) d C M ss ®c N d . We remark that fusion is commutative, so one may swap 
the roles of M and N if desired. 

The proof of the claim (A.23) amounts to showing that any v ® w representing the left-hand side may be 
written as a linear combination of elements of the right-hand side by using the master coproduct formulae 
(A.18). This is demonstrated through the following algorithm which is applied, at each step, to each term v®w 
of the result of the previous step: 


• If v M ss , so v = S n v', with n < —h, then (A.18c), perhaps followed by (A.18b), may be used to replace 
v <8> w by a linear combination of terms of the form v' <8> w' and S m v' ® w, where m > —h. Repeat until none 
of the resulting terms v (8) w have v = S n v' , with n < —h. One may need to take into account relations in M 
to arrive at terms with v G M ss . 

• If w (fi N d , so w = Uw' , where U is a monomial in the chiral modes of weight greater than d, then combine 
A(17) = 0, coming from the left-hand side of (A.23), with (A.18a) and (A.18b) to replace v ® w by a linear 
combination of terms of the form U'v®w' , where the U' are monomials in the chiral modes of weight strictly 
smaller than that of U. Repeat until each of the resulting terms have w G N d , using relations in N. 

• Repeat the above two steps as required. Termination is guaranteed for modules whose weights are bounded 
below over a large class of chiral algebras, the Neveu-Schwarz algebra included, because each step requires 
that the sum of the weights of the factors in each term strictly decreases. 


We remark that these steps are also used when computing fusion products to a given depth d. An explicit 
example illustrating this is detailed in Section 4. 

The first goal in constructing a depth d fusion product (M x N) d is to determine the subspace of M ss <S>c N d 
with which it may be identified. This determination proceeds through the identification of spurious states 
which are actually relations in the tensor product space that are derived from relations in M and N (or even 
in the chiral algebra). Specifically, one combines A(t/) =0, for monomials of weight greater than d, with these 
relations; a spurious state arises if reducing the result to an element of M ss ®c using the above algorithm 
does not yield zero identically. We again refer to Section 4 for examples of this process. Quotienting by the 
spurious states then gives the fusion product to depth d. 

Once the depth d fusion product has been identified, its structure is analysed by computing the action of 
the chiral modes S n with |n| ^ d (all other chiral modes must act as the zero operator). For this, one applies 
the coproduct formulae (A. 18a) and (A. 18b) to a basis element of the depth d fusion product, reducing the 
result to an element of M ss N d using the above algorithm, then to an element of (M x N) d by imposing the 
spurious state relations. By analysing the structures obtained for various (small) values of d, one gets highly 
non-trivial information about the fusion product itself; in favourable cases, the information obtained is sufficient 
to completely identify the fusion product. The example of Section 4 illustrates such a case. 


Appendix B. Staggered Modules and Logarithmic Couplings 

Staggered modules form a particular class of indecomposable modules upon which the Virasoro zero mode 
Lo acts non-semisimply. In a sense, they form the simplest class of such modules and they are responsible 
for the logarithmic structure of most of the best understood logarithmic conformal field theories. The term 
was introduced for Virasoro modules in [53], where some classification results were reported, shortly after the 
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first examples had been exhibited [28]. A full classification of staggered modules over the Virasoro algebra 
was completed in [54]. This notion has been recently developed [48] for other chiral algebras, though only 
a few results have been proven at this level of generality. In this appendix, we discuss the situation for the 
Neveu-Schwarz superconformal algebra and explain the results that we use in the course of identifying the fusion 
products reported in Section 5. 

One particularly important feature of staggered modules is that they are generally not determined, up to 
isomorphism, by the structural diagrams commonly used to depict them. In the simplest examples over the 
Virasoro algebra, one requires an additional parameter ft to fix the isomorphism class completely. This was first 
recognised in [28], but a general invariant definition of ft does not seem to have appeared until [30], where it was 
christened the logarithmic coupling. Other commonly used nomenclature for ft includes “beta-invariant” [54] 
and “indecomposability parameter” [65]. Here, we define logarithmic couplings for the Neveu-Schwarz algebra 
and describe the methods we use to compute these parameters in this paper. 


B.l. Staggered Modules. A staggered module is normally defined to be an extension of higliest-weight mod¬ 
ules upon which Lq acts non-semisimply. Given the results reported here, and in [20,22], we believe that it 
will be necessary to generalise this to extensions of Kac modules. However, most of the rigorous results [53,54] 
about staggered modules were proven for extensions of highest-weight modules over the Virasoro algebra and we 
expect that these results will lift to the Neveu-Schwarz algebra without difficulty. Moreover, the simplest Kac 
modules are also highest-weight modules, so that many of the staggered modules we explicitly analyse will be of 
highest-weight type. We will therefore make use of Virasoro results, lifted to the Neveu-Schwarz algebra, that 
apply to highest-weight type staggered modules, deferring rigorous proofs of these lifts to a future publication. 

Let us therefore define a Neveu-Schwarz staggered module to be an extension of a Kac module by another 
Kac module upon which Lq acts non-semisimply. The short exact sequence characterising the staggered module 
S is then 

0 —» Kr, B -4 K Pl<r —► 0, (B.l) 

where the injection l and surjection tt are Neveu-Schwarz module homomorphisms. We will customarily label 
the staggered module with two sets of indices Sf'j such that the Kac quotient is Kf+k,j+t: and the Kac submodule 
is ICi-k,j-e■ In the above exact sequence, we would therefore write 


c _ < ’5 (p-NA(‘ t - s ) 
3(p+G, 5( ct +W 


(B.2) 


We remark that this labelling need not completely specify the staggered module up to isomorphism. When 
necessary, we shall append the required additional labels in parentheses; for example, (/3). We also mention 

that the staggered modules Sf'j that we have encountered in our fusion computations all have either k = 0 or 
l = 0. However, staggered modules with both k and l non-zero do exist. 

Most of the staggered modules analysed in this paper have the property that both the Kac submodule and 
quotient appearing in (B.l) are highest-weight modules (an example where this is not the case is discussed in 
Section 5.3). In this case, we have the following result. 


Proposition B.l. If the Kac modules K r<s and K p cJ in (B.l) are highest-weight, then their minimal conformal 
weights must satisfy h p ^ a ^ h rjS and h p ^ a — h r>s £ 

If this condition on the minimal conformal weights is not met, then there is no staggered module S making 
(B.l) exact. The requirement that the Kac modules be highest-weight is necessary as the example discussed in 
Section 5.3 shows. In contrast, the following results hold for general staggered modules. 


Proposition B.2. The Jordan blocks of Lq, acting on S, have rank at most 2. 

Proof. Let v belong to a Jordan block for Lq where the (generalised) eigenvalue is h. Then, {Lq — h)v need 
not be zero, but it(Lq — h)v = {Lq — h)irv = 0, since 7 tv belongs to the Kac module /C Pi(T . By exactness, 
{Lq — h)v = lw for some w £ /C riS , hence (Lq — h) 2 v = {Lq — h)tw = l{Lq — h)w = 0. ■ 
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It follows from this proposition, and the definition, that the action of Lq on a staggered module always possesses 
Jordan blocks of rank 2. 

Proposition B.3. Let w,y £ S be elements of a rank 2 Jordan block for Lq satisfying (Lq — h)y = w. If ny 
is annihilated by some U in the universal enveloping algebra of the Neveu-Schwarz algebra, then Uw = 0. In 
particular, if ny € /C PiCT is singular of conformal weight h, then w is singular or zero in i(/C rjS ). 

Proof. We may assume, without loss of generality, that U is homogeneous, meaning that [Lo,J7] = — nU for 
some n € 1. Since nUy = Uny = 0, we have Uy € i(/C r . s ) by exactness, hence Uy is an eigenvector of Lq 
with eigenvalue h — n. Thus, Uw = U(L 0 — h)y = (L 0 — h + n)Uy = 0. The last statement now follows by 
combining this result for U = Lo — h with that for U a positive mode. ■ 

Note that if /C Pj(7 is highest-weight with liighest-weight vector ny of conformal weight h, then w = (Lo — h)y 
cannot be 0 if there are to be any Jordan blocks at all. It follows that, in this case, w is singular. 

B.2. Logarithmic Couplings. In this section, we assume that the Kac submodule and quotient of each 
staggered module is highest-weight. With this restriction, we can follow [48,54] in discussing parametrisations 
of the isomorphism classes of the staggered modules (B.l) with /C r>s and /C PiCr fixed. The parametrisations for 
general staggered modules are beyond the scope of this paper. 

When the highest-weight vectors of the (highest-weight) modules /C rjS and K. PiC , have equal conformal weight, 
h r , s = hp'<r, then this sequence determines the staggered module up to isomorphism (assuming that it exists). 
This follows easily from the fact that such staggered modules are quotients [54, Cor. 4.7] of universal “Verma- 
like” staggered modules. However, the exact sequence typically does not determine the isomorphism class 
uniquely when h p . a > h r _ s . This seems to have been first recognised in the Virasoro examples constructed 
in [28], where an additional parameter /3 was introduced to specify the module structure. The claim that this 
parameter determines the isomorphism class was subsequently demonstrated for a class of Verma-like Virasoro 
staggered modules in [53]. Extending this work to general Virasoro staggered modules required a general 
invariant definition of (3 [30] and was completed in [54]. 

The staggered modules typically encountered in fusion computations have the following structure: Let x 
denote the highest-weight vector of the submodule i(/C r>s ) and let w = Ux denote its singular descendant of 
conformal weight h Pt<7 (unique up to rescaling). The existence of w is guaranteed by Proposition B.3, so we 
may choose y such that ( Lq — h Pt<T )y = w. It follows that U'y must be proportional to x and we define [30] this 
constant to be the logarithmic coupling /3. The choice of y is generally only unique up to adding elements of 
t(fcr,s) (with conformal weight h p ^); however, these are annihilated by as Ux is singular. The logarithmic 
coupling is therefore independent of such choices once we fix a normalisation for the singular vector w. For 
the Neveu-Schwarz algebra, it is not hard to show that the coefficient of the monomial involving only G_i/ 2 is 
non-zero [7], hence we let 

w = (G h _T/ 2 hr ’ s +••■)*• (B.3) 

With the (antilinear) adjoint generated by L* = L_ n and G] = G-j , we see that renormalising ic by a factor 
of a leads to a renormalisation of /? by a factor of |a| 2 . We will often denote the logarithmic coupling of a 
staggered module obtained from fusion calculations, by for convenience. 

Logarithmic couplings are therefore fundamental representation-theoretic quantities. Their physical interest 
in logarithmic conformal field theory lies in the fact that they also appear in the coefficients of certain operator 
product expansions and correlation functions, often accompanied by the factors with logarithmic singularities. 
One can compute the value of any logarithmic coupling by explicitly constructing the staggered module to 
sufficient depth as in the Nahm-Gaberdiel-Kausch fusion algorithm. However, this is, computationally, very 
intensive. An alternative method proposed in [31], but typically limited to staggered modules involving braid 
type highest-weight modules, is to consider staggered modules for which the quotient /C PiCT is replaced by its 
Verrna cover. The desired staggered module may then be realised as a quotient of the Verma-like one if the 
latter possesses a singular vector of the appropriate conformal weight. Checking explicitly for the existence of 
such a singular vector fixes /3 (in the braid type case); see [54] for the proofs. 
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To the best of our knowledge, the most efficient means of computing the logarithmic coupling in a fusion prod¬ 
uct is the field-theoretic, though somewhat heuristic, method outlined in [65] (noting the important clarifications 
described in [103, App. D]). This builds on earlier work [104] addressing the so-called “c -A 0 catastrophe” that 
is reviewed in [105,106]. The computation for the staggered module appearing in (B.l) realises its logarithmic 
coupling /3 as a limit — the parameter t is perturbed away from the desired value t{ 0 ), hence the perturbed 
central charge c and Kac weights h r<s take the form 


t(e) = t{0)+s, c(e) = c( 0 ) — 3M — Je + 


^r,s( £ ) — h-r,s( 0) 


r 2 — 1 S 2 — 1 

8 /( 0) 2 8 


e + 


• (B.4) 


Let x(e) denote a highest-weight vector of central charge c(e) and conformal weight h r ^ s {e). We define U such 
that Ux( 0) is the singular descendant of conformal weight h UyV ( 0), normalised as in (B.3), then let w(e) = Ux(e). 
Note that U does not depend on e. The logarithmic coupling is then given by [65] 

8/(0 ) 2 (x(e),U^Ux(e)) 

inn 


/3 = — lim 


{w{e),w(£)} 


(B.5) 


■ mu - - — ---;— -—--— ----—— — - 7T inn 1 

£— h UtV (e) — h r ^{e) — (h U:V ( 0) — /i r , s ( 0)) u 2 — r 2 — ( v 2 — s 2 )t(0) e_! '° 

where x(e) is normalised so that (x(e), x(e)) = 1. We provide an example illustrating the use of this formula in 
Section 4. 


Appendix C. Further explicit fusion computations 

In this appendix, we list a selection of the fusion products that we have computed by combining the informa¬ 
tion provided by the Verlinde formula with explicit Nahm-Gaberdiel-Kausch calculations and staggered module 

theory. As it can be difficult to verify explicitly that a sum is direct (the depths required can be very large), we 

? 

introduce a symbol “©” for these fusion rules to indicate that the sum indicated may, or may not, be direct. We 
do not indicate parity — see Section 5 for this information. All logarithmic couplings have been independently 
verified using (B.5). 


(p,p') = ( 1,3) (c= -§). 


1 C 1,3 

X 

1 C 1,3 = 

= <Si 0 3 2 ffi/Ci, 3 , 

/C 3 ,i 

X 

/C 3 ,i = 

/Ci,i 

© 

1 C 3 , 

/Cl, 3 

X 

/Ci ,5 = 

= /Cl ,3 ©<’6 (-2), 

/C 3 ,1 

X 

/c.5,1 = 

/C 3 ,1 

© 

/c 5 , 

/Cl, 3 

X 

/Cl, 7 = 

= 5 1 % 1 (- 2 )©/Ci, 9 , 

/C 3 ,1 

X 

/c 7,1 = 

/c.5,1 

© 

1C7, 

/Cl ,3 

X 

/Cl ,9 = 

= 5 ^(- 64 / 9 )©/Ci, g , 

/C2,2 

X 

/C 2,2 = 

/Ci,i 

© 

/Cl, 

/Cl ,3 X /Cl, 11 - 

= 5 ?;i 1 2 (- 64 )©/Ci, 9 , 

/C2,2 

X 

/C 4 ,2 = 

/C 3 ,1 

© 

/c 3 , 




/Ci, 5 x /Ci,5 = 7 ^.°[ 2 

©/c 

1,3 

© /Cl,7 

©/C; 

1,9 



/C5.1, 

^7,li 

^9,li 

^ 3,1 © ^ 3 , 3 , 

/C. 5,1 © /C5.3, 


(C.l) 


(p,p') = ( 1,5) (c=-§). 

JCl ,3 X JCl ,3 = /Ci,i © /Cl,3 © /Cl,5, 

/Cl,3 X /Cl ,5 = Cq,’ 2 © /Cl,5, 

/Cl ,3 x /Cl,7 = /Ci,5 ® /Ci,7 ® /Ci,9, 

1C 1,3 X K 1,15 = Si°; 2 5 (- 1152/25) © /Ci,i 6 


/Ci, 5 x /Ci, 5 = S 0 ^ © ^ © K i, 5) 

/Ci,5 x /Ci, 9 = 5 1 °; 1 3 0 (16/25) © /Ci,5 © ^(-d) 
/Ci,3 x /Ci,9 = /Ci,7 © »S 3 ’j L g(— 4 ), 

/Ci, 3 x /C 1,19 = 336) © /Ci, 17 . 


(C.2) 
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(P,P') = (2,4) (c = 0). 


IC 1,3 

X 

/Ci, 3 = 

II 

© 

O- 

1), 

/C 3> 1 

X 

/C 3 ,1 

1C 1,3 

X 

/Ci, 5 = 

c 0,i 
“ D l,4 

(- 

1 )©/C 

1 , 7 , 

/C 3 ,1 

X 

/C5,l 

/Cl,3 

X 

ic 1,7 = 

= fcl,5 

© 

O- 

15), 

/C 3) 1 

X 

IC 2.2 

/Cl ,3 X /Cl ,11 - 

= ^T,9 

© 

S°iU- 

-3780), 

/C 3 ,1 

X 

/C 2 ,4 

/C 1,3 

X 

IC-2,2 = 

= ^2,2 

© 

!C 2 , 4 , 


ic 1,5 

X 

/Cl,5 

IC 1,3 

X 

/C 2 , 4 = 

c 0,2 
“ °2,4 

© 

IC 2 , 4 , 


/Cl, 3 

X 

/C2,6 


IC 2,2 X ^-2,2 = «5>2,’l (3/8) © 


= 5 2 1 ; 1 °(3/8)0/C 5> i, 

= 5 4 1 ; 1 °(14175/32)©/C 7> i 
= /C;2,2 © ^4,2, 

= /C2,4 © ^4,4, 

= 5 1 °; 4 3 (-1/4)©5 1 °; 4 1 (-1)©/C 

= /C2,4 © ^2,6 © ^-2,8, 


1 , 9 , 


(C.3) 


(2,8) (c 

= — 

21 

4 

)• 









/Ci,3 x 

/Cl ,3 

= 

IC 1,1 i 

B /Ci, 3 © /Ci, 5, 

ic 1,3 

X 

/C2 ,2 — 

/C 2 ,2 

© 

/C 2 ,4 


/Ci, 3 x 

/Cl ,5 

= 

IC 1,3 € 

B /Ci, 5 © /Ci, 7, 

/Cl ,3 

X 

/C 2 ,4 = 

^2,2 

© 

/C 2 ,4 

© /C2,6> 

/Ci, 3 x 

/Cl ,7 

= 

0 - 

-3) ©/Ci, 5 , 

/Cl ,3 

X 

1 C 2.3 = 

/C 2 , 4 

© 

/C2,6 

© /C 2,8, 

/Ci, 3 x , 

*Cl ,15 

= 

/Ci , 13 

©<18 (-165), 

/Cl ,3 

X 

IC 2,8 = 

c 0,2 

°2,8 

© 

/C2,8i 


X 

CO 

*Cl ,23 

= 

/Cl, 21 

© t S 1 °’ 2 1 4 (-163020), 

/C 2,2 

X 

/C 2> 4 = 

c i,0 

° 2,3 

(- 

105/256) © 5 2 ’5 (—15/64), 





IC 2,2 x /C3,3 = IC 

2,2 © /C2 

,4 © /c 4 , 

2 © IC 4, 

4 - 



5) (c 

= 

-). 
10 > 

/Ci ,3 

X 

ic 1,3 

/Ci ,3 

X 

/Ci, 5 

/Cl ,3 

X 

/Ci ,7 

/Cl ,3 

X 

/Ci, 9 

/Cl ,3 

X 

/C 2,2 

ic 1,3 

X 

/C 2 ,4 

/Cl ,3 

X 

/C 2,6 

/Cl ,3 

X 

/C 3,3 

/Cl ,3 

X 

/C 3,5 


£1,1 © /Ci ,3 © /Ci,5, 

5 i °’.5 (- 256 / 675 ) ©/Ci >5 , 
/Ci,5 © /Ci,7 © /Ci,9, 

/Ci,7 © ‘S'iJio(— 11264 / 9 ), 

K. 2,2 © /C 2 , 4 , 

/C 2 ,2 © <52,5(-2/3), 

<5 (-2/3) © /C 2 ,8, 

/C34 © /C 3,3 © /C 3 , 5 , 

S 3 ,5 © /C 3 ,5, 


/C 3 ,1 

X 

/C 2,2 = 

c i,o 
“ d> 3 2 

(64/125), 

/C 3 ,i 

X 

/C 4> 2 = 

c i,0 
“ d> 3 2 

(64/125) ©/C 6 , 2 , 

/C 3 ,i 

X 

IC 3.3 - 

c 2,0 
“ °3,3 

(256/1125) ®/C 3 , 3 , 

/C 3 ,i 

X 

ic 2 ,4 - 

c i,0 
“ °3,4 

(2/5), 

/C 3 ,1 

X 

IC 3,5 = 

c 2,0 
“ °3,5 

© /C3,5, 

/C2,2 

X 

IC 2.2 - 

= /c M 

© /Cl ,3 © /C3,1 © /C3,3 

/C2,2 

X 

ic 2 ,4 - 

= /c 1>3 

© /Cl ,5 © /C3,3 © /C3,j 

/C2,2 

X 

IC 3,3 = 

c i,0 
“ ^3,2 

(64/125) © C>3 4 (2/5) 

/C2,2 

X 

/C 4 ,2 = 

= JCs, 1 

© IC 3,3 © 1C. 3,1 © /C 5,3 


(C.5) 
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